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Preface 

In this study, the authors have attempted to lay the groundwork need¬ 
ed to determine the feasibility of a trajectory-match technique for inter¬ 
cepting intercontinental ballistic missiles. In Chapter II, the ballistic 
trajectory equations end rocket equations of motion are developed in de¬ 
tail, since the authors feel that many readers, who may be interested in 
I.Q3M defense, may not be completely familiar with these concepts. 

We found the boost-phase analysis in Chapter IV to be the most 
challenging aspect of the problem. As far as we know, we have developed 
a new technique for obtaining good approximations to the thrust and 
weight of boost-rockets required for boosting a given payload onto a 
specified trajectory. 

Figures are used liberally throughout the report so that the reader 
can more easily follow the geometry involved in some of the more compli¬ 
cated mathematical derivations. There is also a list of definitions 
included on page vl to explain some terms which are peculiar to this 
particular study. The computer programs used are included in Appendix A, 
so that other investigators nay reproduce any of the results or expand 
this study, if they so desire. 

Ve wish to thank Captain H.G. Pringle for suggesting this thesis 
topic, and for his many helpful ideas throughout the investigation. 

The authors are also indebted to their wives and children, for 
their patience, understanding, and many sacrifices which have greatly 
contributed to the completion of this work. 

H. A. C. 

C.A.M. 
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Definitions 


Interceptor Vehicle The final payload consisting of the 
vehicle structure, and equipment necessary for detection and 
destruction of the enemy warhead. 

Velocity-Match Rocket An assembly consisting of the interceptor 
vehicle, rocket engine, fuel tanks, eto., needed to accelerate the 
interceptor vehicle to its final velocity. The velocity match 
rocket is the payload for the boost-phase rooket. 

Interceptor System The structure consisting of the boost-phase 
rocket and the velocity-match rocket. 

Excess Time The difference between the flight-times of the 
• * 

enemy warhead and the interceptor vehicle from their launch sites 

4 

to the intercept point. 

vS* - 

Sub-Bocket An assembly consisting of the payload and the 

.V 

^'rocket stages not yet jettisoned. 
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cross sectional area of rocket stage 

aiming point 

speed of sound 

height above the earth 

coefficient of drag 

effective mass ratio of "boost-phase sub-rocket 
effective mass ratio of velocity-match rocket 
drag force 

total energy per unit mass 
energy ratio 

eccentricity of elliptical trajectory 
functional notation 
acceleration of gravity 

2 

proportionality constant (lb mass per lb force)(ft per sec ) 

angular momentum per unit mass 

specific impulse 

intercept point 

kinetic energy per unit mass 

number of stages 

launch site 

initial 3ub-rocket mass (lb^), numerically equal to sea level 

weight 

mach number 


vii 
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m stage mass (lb ) 

m 

nip interceptor vehicle mass 

nip interceptor rooket mass 

n structural factor 

r length of radius vector from center of the earth 

p 

I radius at intercept point 

p 

la radius at aiming point on corrected intercept trajectory 
r^ radius to interceptor trajectory 

r^ radius to enemy warhead trajectory 

S distance 

S' distance travelled by one-stage rocket, starting from rest in 

field-free, airless space 
thrust force of the jth stage 

T* flight time from apogee along ballistic trajectory 

T x excess time (See list of definitions) 

Th. ratio of thrust of jth stage to weight of interceptor vehicle 

J 

t time 

u rocket exhaust velocity 

V potential energy per unit mass 

v velocity 

range over earth's surface 
x horizontal distance 

y vertical distance 

a angle between thrust vector and velocity vector 

p thrust to initial weight ratio of boost-phase sub-rocket 

thrust to initial weight ratio of velocity match rocket 


viii 
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Y 

r 

6 

€ 


V 

v 


+ 1 


4 

<t> 

<*> 

4 > 


Lla 



mo 

moa 


P 


f 

V'u 

X 


flight path angle (angle between velocity vector and local 
horizontal) 

flight-path angle at termination of tip-over maneuver 

angle between aiming point and intercept point 

angular difference between launch site and burnout point in 

rectangular coordinates 

gravitational constant of the earth 

marimim acceleration in sea-level w g’s" 

burnout angle in rectangular coordinates 

angle between thrust vector and local horizontal 

position angle between radius vector and major axis 

burnout range; angle of corrected interceptor trajectory 

angle from launch site to intercept point for theoretical 

intercept trajectory 

angle from launch site to intercept point for corrected 
intercept trajectory 

interceptor launch site position relative to warhead apogee 

semi-range angle of theoretical intercept trajectory 

semi-range angle of corrected intercept trajectory 

atmospheric density 

stage mass ratio (boost phase) 

overall mass ratio (boost phase) 

overall mass ratio (velocity match rocket) 

thrust angle above horizon in rectangular coordinates 
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corrected values 

final value 

burnout conditions o 
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vector notation 
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Abstraot 

In this thesis, a mid-course, trajectory-match intercept of 
intercontinental ballistic missiles is analyzed. The intercept 
is accomplished by using a multi-stage rocket to boost an inter¬ 
ceptor vehicle onto a ballistic trajectory in the same direction 
of flight as the warhead trajectory, so that the two trajectories 
are almost tangent at a selected intercept point. Near the 
intercept point, a final rocket is ignited to accelersti the 
vehicle and match the warhead's position and velocity. 

A simplified analysis is performed for the case of a co-planar 
intercept over a spherical, non-rotating earth. Intercept 
trajectories are determined from the trajectory geometry and the 
rocket equations of motion. An approximate method is derived to 
estimate the thrusts and rocket weights required to boost an inter¬ 
ceptor vehicle onto these trajectories. Equations are programmed 
for an IBM 1620 digital computer to obtain estimates of thrusts, 
weights, maximum accelerations, and the difference in the times of 
f flight of the warhead and interceptor. Various combinations of 
intercept points and interceptor launch sites are used to determine 
the first stage thrust required per pound of interceptor vehicle for in¬ 
tercept of low angle, minimum energy, and high angle, 6000 and 4000 
nautical mile range warhead trajectories. The effects of specific 
impulse, structural factor, and thrust to initial weight ratio upon 
overall thrust, time, and acceleration are also determined. It is 
found that there is an optimum value of thrust to weight ratio which 

xiv 


GAW Mech 61-2 


results in minimum overall thruBt requirements. 

A numerical solution to the equations of motion is performed on 
the computer for the "boost phase. For the case considered it is 
found that: the actual "burnout velocity is within three percent 
less than desired and thrust is within eleven percent of that 
required, thereby proving the validity of the original "boost phase 
approximations. 

The computer results are presented graphically, and indicate 
that this intercept technique is theoretically possible for either 
a manned or unmanned interceptor system. However, a much broader 
investigation would be required to definitely determine its useful¬ 
ness for defense aga?^st ICBM's. 


xv 
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AN ANALYSIS 
OF A 

TRAJECTORY AND VELOCITY MATCH TECHNIQUE 
FOR INTERCEPTING 

INTERCONTINENTAL BALLISTIC MISSILES 

I. Introduction 

General Considerations 

The high speed and short time of flight of the intercontinental 
ballistic missile have made the problem of finding an effective ICBM 
defense system extremely complex. To successfully intercept an ICBM, 
an interceptor vehicle must be positioned, at or near the same i'Ant in 
space as the ICBM, at the same time. The interceptor must also be 
able to distinguish the warhead from decoys and destroy it. 

One system which is presently in the development and test stage i - 
limited to intercepting enemy missiles during the re-entry phase of their 
trajectories. For this type of intercept, the interceptor vehicle 
approaches the warhead praotioally "head on", resulting in extremely high 
relative velocities. This places very stringent guidance and accuracy 
requirements upon the system. Furthermore, this system is limited to 
pomt-aeieaae of relatively small, strategic target areas. Since tin. i4» 
practically a "last ditch" defense, the question has arisen: Can a 
successful intercept be made during the boost phase of the warhead 
trajectory, while the rochet motors are still burning, or during the 
mid-course, free-flight phase, prior to re-entry? 


1 
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It would be difficult to intercept a missile during the boost phase 
because of the very short warning time available, and the usual lack of 
friendly defense bases in the vicinity of enemy launch sites. Thus, 
serious disadvantages also exist for this type of intercept. 

Now let us consider the possibility of a mid-course intercept, c<e£. 
a method of reducing the relative velocity between the interceptor and 
the warhead. This can be accomplished by: launching an interceptor 
vehicle on a ballistic trajectory in the same direction of flight as 
the warhead; making this intercept trajectory nearly tangent to the 
warhead trajectory; and accelerating the interceptor vehicle to match 
the warhead's position and velocity by igniting a final booster rocket 
r.s:-.,r the intercept point. This method offsets the disadvant gcs 1 ihc- rt 
in the re-entry and boost phase intercepts. Available warning time would 
not be as great as for the re-entry intercept, but terminal guidance re¬ 
quirement a should not be nearly as critical due to the lor relative - r elo~ 

Jxnear the intercept point. Thus, the mid-course, trajectory-match 
intercept appears to be a promising solution to the problem ul iCBh 5; , - ~.-nse. 

Statement of the Problem 

The purpose of this thesis is to analyze a mid-course, trajectory- 
match 5.ntercept of an I CBM warhead. This analysis consist* of determining 
inhere action geometry, the required thrust and weight of -V . 
esptor system, and the difference between the flight times of the ICBM 
warhead and the interceptor, from launch site to point of intercept. 

Since for this type of intercept there is a low relative velocity 
at the intercept point, there are possible advantages to having a man 
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in the interceptor vehicle. These include man'a ability to mate decieit 
and evaluate data from abort range, decoy discrimination equipment. 
Therefore, the analysis includes a consideration of the maximum 
acceleration imposed upon this system, and sin effort is made to keep 
acceleration within human tolerances. 

Analysis of the Problem 

The analysis of an ICBM intercept over the non-spherical, rotating 
earth, by an interceptor system of unknown performance would be extremely 
complex. However, this analysis can he greatly simplified by assuming 
that the earth is spherical and non-rotating, and that the warhead and 
interceptor trajectories are co-planar. Although these assumptions are 
unrealistic, the basic concepts and techniques, when established fox tluu 
-necial case, could be extended to the general intercept, .wii -» 
^resent investigation show promise- 

Interoept Analysis . The first step is the detersa? na. io-> ' p th'i 
geometry of the intercept for the tangent trajectory. The equations 
are developed to allow variable inputs of warhead trajeotory, intercept 
point and interceptor launch site. 

Next, the dynamics of the velocity-match phase are considered. It 
is shown that the equations of motion of a rocket in field free space can 
v e considered to be a perturbation to the theoretical free-1‘3 xght 
trajectory, if gravity is assumed constant during the velocity-match 
phase. The required free-flight trajectory is determined by correcting 
the theoretical trajectory for the effect of this perturbation. 

Boost Analysis . With the required intercept trajectories 


3 



GAW Mech 61-2 


established., the next problem is to determine the booster requirements 
to place the interceptor vehicle on this trajectory. These require- 
ments must be based upon the weight of the velocity-match rocket which 
in turn depends upon the weight of the interceptor vehicle. Since this 
is an initial investigation, the weight of the interceptor vehicle is 
unknown. A vehicle weight could be assumed, but to permit greater 
flexibility, non-dimensional weight ratios are used to determine rocket 
weight and thrust requirements per pound of interceptor vehicle. The 
boost phase analysis is further complicated by the large number of 
variable rocket parameters. The thrust, weight, specific impulse and 
structural factor of each stage, as well as the number of stages must be 
established to continue the analysis. Furthermore, the exact sire \~J 
shape of the interceptor system and the velocity profile during bo^t 
all unknown. Therefore aerodynamic drag during the boos+ phase oar r~.+ 
be exactly predicted. These complications make a closed-form solv.ti -- 
to the boost-phase equations of motion unfeasable. 

The following assumptions are made to obtain an approximate solution 
to these equations: 

(l) The stiuctual factor and specific impulse are the same for each 
stage. 

{2) Drag is negligible compared to thrust. 

( 3 ) Gravity is constant in magnitude and direction. 

( 4 ) Vertical launch and tip-over maneuver can be neglected. 

(5) An inertial, rectangular coordinate system can be used. 

(6) Thrust for each stage is constant. 
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(7) The thrust direction is the seme (constant direction relative 

to inertial spaoe) for each stage. 

A further simplification is made by making the thrust to initial weight 
ratio the same for each sub-rocket. 

With the use of the above assumptions and simplifications, a ev.i, 
of equations are obtained which give the thrust and weight ratios 
required to achieve a desired velocity and direction in spaoe. These 
equations are then combined with the interoeptor trajectory equations 
in a program designed for use on the IBM 1620 digital computer. By 
means of this program the thrusts and weights of the booster rockets, 
the flight times from interoeptor launch site to intercept point, and 
maximum accelerations are obtained. 

A numerical integration is then performed on the 1620 computer 
,w check the validity of the boost-phase approximations. For inc 
•ypical intercept situation, it is found that the burnout velocity oi 
tue booster rockets differed from that desired by approximately three 
per cent. 

Effect of Variables. The large number of variables require that 
this investigation be limited to a few of the many possible enemy 
warhead trajectories. 

Two typical ICBM ranges (6000 and 4000 nautical miles) are selected 
Heaver, this must be further limited since for any particular varhoa,' 
range there are an infinite number of possible trajectories, characterised 
by the height of their apogees. In an attempt to adequately consider 
possible intercept conditions, a high, medium (minimum energy), and low 
altitude trajectory is considered for eaoh range. 
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The effect of etructual factor, specific impulse, and. thrust to 
weight ratio on thrust requirements, excess time, and maximum 
acceleration is investigated for a few typical intercepts. Finally, 
reasonable values of specific impulse and structural factor, and ft- 
optimum value of thrust to weight ratio are used to investigate the 
effect of varying the intercept point and interceptor launch site on 
thrust and excess time. 
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II. Basic Theory 

The Free-Flight Missile Trajectory 

Basic Equations . The motion of a Ballistic missile in a vacuut. 
under the influence of a central force field such as the earth* 0 yravi-- 
tational field, can he analyzed hy considering the missile to he a 
point mass in a polar, inertial coordinate system (Ref 13:2-20). 



Let i he a unit vector in the r direction and j a unit vector per¬ 
pendicular to r as shown in Pig. 2-1. The derivatives of the unit 
vectors with respect to 4> are 


dt 
d t 



and 



Then the derivatives with respect to time are 


di _ di dt 
dt d4> ait 



and 


— - d>T 

dt dt dt 


7 




GAW Meoh 61-2 


Now the position vector of the missile is 


Its velocity is 


r= re 


r = r i + r 


r - ri + r<P j 


and its acceleration is 


or r =(r - r<p Z j L -h{2.rcp -h rep) j 

How ^r(r 2 <p) ~ 2rr<p -h r z cp = r(Z r<j> + r<p) 


~r~ Jt(r 2 <p) = 2rj> + rp> 


( 2 - 1 } 


( 2 - 2 ) 


(2-3) 


Substituting Eq (2-3) into Eq (2-2), 


r-(r-ri‘)i +[yr 3r(r‘i)]j 


Now, by applying Newton's Second Law of Motion, P ■ ma 


where 


£7 m — 

F=-~pr- l 


m is the mass of the missile 


p is the earth's gravitational constant 


we get 


Mm ~ 


L ~ m[(r- r <p z )i -h(yr ^t[r z <p]),'jj 


Then by equating the j components of Eq (2-5) 

I d / z ; \ - 


~’dT( r ^*) = ° 


(2-5) 


Now, since r is not infinite. 
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lit ( r& 4 > ) = 0 


or r a cfi> — "a constant" 

But r <b is the polar coordinate expression for the angular momentum ; 
unit mass. Therefore, the angular momentum of the missile is a constant 
during the free-flight phase. Let us call the angular momentum par ~-r 1 
mass, h. Then 

h=r z 4> ( 2 - 6 ) 

Now, equating the I components of Eq (2-5)> 

- ~ r - r<j>* (2-7) 


In order to solve Eq (2-7), let us define a dummy variable, 


Then, in terms of p, Eq (2-6) becomes 


h=r*4‘ = jz 


or 


Also 


4> = ph 


But 


d<f> 

_ L 

d4>\ P J ~ P* d<f> 

Combining Eqs (2-9) and (2-10), 


Then 




d<p 



(2-B) 


(2-9) 


( 2 - 10 ) 


/ - \ 
V ^ * ' 


9 




OAW Meoh 61-2 


( 2 — 1 ?) 


or 



d Z p 

d<t> 9 


Substituting Eqs (2-8 and (2-12) into Eq (2-7), with r ■ ^ , we get 

2 p 

Dividing Eq (2-13) by -pa gives 

(~M/h*) = -t p ( 2 - 14 , 

Equation (2-14) is a non-homogeneous, second order, linear differential 
equation which can be solved by the method of undetermined coefficients 
(Ref. 14*36-41). 

The solution of this equation is 


p = jp: -h A COS (4>~ <t>') (2-15) 

where A and <b ' 416 constants which must be determined from the boundc. 
conditions. 

The evaluation of these constants depends upon the feet that thr 
.otal energy of the missile must be a constant for free-flight in a 
central force field. Thus, 

E - K + V - "a constant" 

where: E » total energy of the missile per unit mass 

K - kinetio energy of the missile per unit mass 
V - potential energy of the missile per unit mass 
By defining potential energy, V, as zero when r is infinite we can v*.y 


that 


V= - 


r 


Also, kinetio energy is defined as 



or 


K=2~(r z -t- r*<p z ) 


( 2 - 16 ) 

( 2 - 17 ) 
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Thus 


s') M 


E-±(t-' + rV)- 


Equation (2-18) can be rewritten in terms of p by the use of Eqs ( 2 - 3 ' 


and (2-8). 


£ ~ ~k[h z (ijf) +- (~p)(p z hf] ~JJ P 


h * rsdp\i 


2 Hd 


V + P j ~^P 


(2-...V) 


Differentiating Eq (2-15) gives 


^if = - A Sin (4>- c P') 


( 2 - 20 ) 


By using Eqs (2-15) and (2-20), p is eliminated from Eq (2-19» giving 
an expression for E in terms of A and h. 


h 2 A a 

£ = a 2h z 

Rearranging Eq (2-19a) and solving for A, 


Substituting Eq (2-2l) into Eq (2-15)» 


( 2 - 19 a) 




+ * -V ,, 2Eh z 




h‘~ h 


p ~j?[' cos c*-^‘)] 


( 2 - 22 ) 


J_ = h%U _ 

P cos (^~^) 


(2-v' 


Now if we choose our coordinate system so that r is a maximum when 
<p • 0, then 4 >/ m 0 and Eq (2-23) becomes 
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(2-24 J 


h*/M 

I - V I + -^fr cos 4 


This equation is similar to the equation of an ellipse in 
polar coordinates, which is 


t-ecosf ( 2 - 25 ) 

where 1 o the semi-lattice rectum of the ellipse 

e - the eccentricity of the ellipse 
Therefore we can state that the equation of the trajectory of a ballis¬ 
tic missile in free-flight over a spherical non-rotating earth is an 
ellipse with the center of the earth at one focus. The eccentricity of 
the ellipse is defined by the equation 



and the semi-lattice rectum by 





(2-27) 


Thus the trajectory of a ballistic missile is completely defined if 
the total energy and angular momentum of the missile are known, and is 
given by the equation 

hVM 

/-s cos 4> 

Additional Trajectory Parameters. Figure 2-2 shows the elliptical 
trajectory plus several other parameters which are useful for analyzing 
ballistic trajectories. These parameters are defined as follows: 

Y - the angle between the flight path and the local horizontal 

v - theoretical launch angle at the earth's surface 

1 0 
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v q - theoretical launch velocity at the earth's surface 

X - range over the earth's surface from theoretical launch 

© 

site to theoretical impact point 

r - radius of the earth 
o 

4> - semi-range angle of the trajectory 

T 0 

4> - angle from major axis of the ellipse to the radius 

vector (measured positive clockwise) 

The angle of the tangent to the trajectory, can k© determined 
from the ellipse equation and Fig. 2-3. 
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tar) r - c/r = 1 dr 
ion r - - r d< t> 

(2-2;? 

Differentiating Eq (2-25), 


dr -1 e sin 4> 

d<t> (l-e cos <£)* 


Then 


tan r- ,' eW , 

1 ~ <2 COS <p 

(2-31) 


Prom Eq (2-31), the theoretical launch angle is defined "by the equation 



Tiae of Flight . One very important parameter of the trajectory 
remains to he determined. This is the time of flight during the free- 
flight trajectory (Ref. 13»17-19)* An inspection of Pig. 2-3 shows that 
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the incremental area between two radial lines, d<£ apart, is 


dA = ~^(rd4>)(r) « r z d4> 


Equation (2-36) can he written as 


^ . I id* , c . 
d t ~ 2. r dt ~ 2 r ? 


I 


hut 


= /? 


(2-3.6) 


(2-6) 


and so 


dA 
dt ~ 2 


= ~5~ h 


(constant) 


(2-37) 


For one complete orhit or revolution, 

A=~^hT (2-38) 

where A • the total area of the ellipse 

T - the period of revolution 

Since the time rate of change of the area of a sector of the ellipse it» 
a constant, 

T*_ A* 

T A (2-39) 


idiere T* is the time required to travel along the arc of a sector of 
area A*. 


Then 


T A It 


But, hy rearranging Eq (2-38), 


A * 

T _ 
A 


4 r , i 


2' 


Thus 


T*=fA* 


(2-41) 


How, since A* is the area of any sector of the trajectory, let us 
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evaluate A* for a sector from the apogee radial line (where 
to the general angle , 4> , downrange from apogee. 

From Eq (2-36) . 

dA *-j- r z d<f> 

Squaring Eq (2-25) gives 

Z a _ 

r ~ (/- e cos <f)* ^ A (l~e cos <p) z 
Combining Eqs (2-36) and (2-42) 


Then 


dA = 


h 4 d4> 

2^u 2 (l- e cos <£)* 


A > h* f * d* 

J o 0-e ccs 4>) z 


- o) 


( 2 - 42 ) 


(2-43) 


(2-44) 


This integral may be evaluated by the use of a table of integrals 
(Ref. 8t41)42). The result is 



/ h 4 ) f I \f ~ e sin4> 
2yu z /\l- eVjf /-e co5 



, ~iMl-*' tan 
ian (—— 




Substituting Eq (2-45) into Eq (2-41) we obtain the equation for the 
time of flight from apogee to any angle, » downrange from apogee. 


T 


{Mdpfc 


e st n 4> 


e cos 4 




-c* tan 
/-e 



(2-46) 


Since the trajectory is symmetric about the apogee radius (or nol 
axis)) Eq (2-46) can also be used to determine the time to apogee from 
any angle up-range by substituting the absolute value of the position 
angle in the formula. 

The time of flight between any two points on the trajectory is 
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X*- T*M - TU) 

(2-47) 

where <£,> <f>^ and the correct signs are used for the position angle, 
Minimum Energy Trajectory. (Ref. 13*1924) Previously? we have 
shown that the trajectory of a ballistic missile is completely defined 
by the total energy, E, and the angular momentum, h. Now we will 
investigate the effect of each of these parameters on the trajectory. 

First, let us note from Eq (2—48) that sb the total energy increases, 
the theoretical launch velocity, v, increases. 



-L V Z — 
a v ° 


r e 


(2-48) 


But 


h= r 0 v 0 cos r 0 


(2-33) 


Thus for a constant launch angle, the angular momentum also increases 
with an increase in total energy. 

Rearranging Eq (2-28) gives an equation which relates the angular 
momentum and eccentricity to the semi-range angle 


COS <t> 0 



Substituting Eq (2-26) into Eq (2-49) gives 


(2-49) 



(2-50) 


Inspection of Eq (2-50) ehows that as E increases with h held 
constant, cos <fi decreases or <fi Q , and therefore range, increases 
(See Fig. 2-4a). 

Examination of Eqs (2-33, 48, and 50)shows that for a constant 
launch angle, an increase in the total energy increases both the range 
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and the angular momentum (Pig. 2-4h). 

Now let us determine the effect of angular momentum on the total 
energy required for a fixed range. Equation (2-50) can he written ?? 



If there is a minimum energy for a given range, we can determine the 
value that h must have for minimum energy hy differentiating Eq (2-5l) 
with respect to h, equating the result to zero, and solving for h. 
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Finally, for a given range there is an almost unllmit-.,i 
trajectories possible depending on the values of E and h. For any 
particular value of E which is greater than E^^, there are two possible 
trajectories for the same range. These two trajectories converge to t>? 
minimum energy trajectory as the total energy approaches E ^ (See Fig. 
2-4d). 

Now that we have established the existence of a minimum energy 
trajectory, let us derive a few minimum energy relationships. Substi¬ 
tution of Eqs (2-54) and (2-55) into Eq (2-26) gives the eccentricity 
of the minimum energy trajectory as 

C os <£,_ 

( 2 - 56 ) 


& — 
'— M P" 


ME - I -h sin 4> 

The optimum launch angle is determined by substitution of Eq (2-56) 
-to Eq (2-32). 


tan x 


_ cos 4> o 


5 m.£ /-f sin 4* 

1 o 

Notice for the minimum energy case, 


(2-5?) 


e. tan X MJE 


Energy Ratio Concept . The energy ratio (Ref. 6: 134) is a 
measure of the energy available for the trajectory, and i3 defined as, 

r a V B 


ER = 




If the range, or <p> and the energy ratio are specified, the launch 
angle needed to achieve this range is determined. The equation is 

tan r s = r/rlfx (2-58) 
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This equation is used to specify the high and low angle trejec:t' :’ic ;■ 
which are investigated in the latter portion of this report* 


Rocket Equations of Motion 

Velocity Equations . The motion of a rocket can "best be describe! 
by considering the case of a rocket in field-free space (no external 
forces). Let us assume that we have a rocket in this environment, 
traveling horizontally with velocity, v, and-mass, M, at time, t. 

(Fig. 2-5). 

E > ~* v 

Rocket at time * t 

u-v M-mt ^ —► v+dv 

Rocket at time, t + dt* 

Fi $ . 2-JT 

Diagram -for Rocket Equation Derivations 


Let t. 

o 

M 



s 

0 


% 

m 


u 

Since no 
momentum must 


= the total burning time of the rocket 

■ the mass of the rocket before ignition 

= the mass of the rocket at burnout 

■ the velocity of the rocket when t * 0 

= the velocity of the rocket when t * t^ 

= the displacement (position) of the rocket at x *» ■.’ 

a the displacement of the rocket at t = t^ 

* mass rate of flow (burning rate) of fuel, a constant 

* the exhaust gas velocity relative to the rocket 
external forces are acting on the rocket of Fig. 2-5, the 
be a constant for the rocket system (including exhaust 
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(2-69) 
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s-s.i-v.t + ujjnpi^dt 

(2-?^; 

u M T Mm 

Hence ,S- s o ~h V 0 t +■ ut . In —rr 

O o m A7 

(2-75) 

For the burnout conditions, Eq (2-75) becomes 


S b - S a + V 0 t k -b S' 

(2-76) 

wh9r0 , , uM f 7 M. 

(2 7 *' 

or 5= u[t t --%lnD] 

(2-78) 


Equations (2-70, 76, and 73) are the "basic rocket equations which 
used in the analysis of the problem. They can be modified as necessary 
to account for the effects of gravity and atmospheric drag. 
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Specific Impulse of Rocket Fuels . The velocity and distance 
equations of rocket motion are both functions of the velocity of the roo 
exhaust gases. Usually, the exhaust velocity, which is a characteristi 
of a given fuel, is not stated explicitly. Instead, the specific 
impulse of the fuel is used to define the exhaust velocity by the 
equation (Ref. 10 s 20) 

u » I s g c (2-79) 

where I , the specific impulse, is the pounds of thrust per pound of 
s 

mass flow per second, and g is a proportionality constant numerically 

c 

equal to the sear-level acceleration of gravity. 

The derivations which follow in Chapters III and IV have 
exhaust velocity, u, as a variable. However, the computer progrsssp up* 
to solve the equations have I # , rather than u, as an input variable. 
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III. Determination of Intercept Trajectory 


Theoretical Intercept Geometry 

As a first step in determining the interceptor trajectory let us 
compare two possible cases as shown in Pigs. 3-1 and 3-2. In the 
"apogee" intercept, the point of intersection between the warhead 
and intercept trajectories occurs when the interceptor is at the 
apogee of its trajectory. On the other hand, for the "tangent" inter¬ 
cept, the intersection is at the common tangent to both the warhead 
and intercept trajectories, and this point is "down range" from the 
interceptor's apogee. By comparing Pigs. 3-1 and 3-2 it is seen that 
for the same intercept point and interceptor launch site, the tangent 
intercept will result in a smaller vector velocity difference between 
the warhead said the interceptor. This is an important fact, since 
the weight of the velocity-match rocket increases exponentially 
with the magnitude of this velocity difference. Therefore this 
investigation is limited to the tangent intercept technique. 

The most direct approach would be to fix the desired intercept 
point and interceptor launch site, and then determine the intercept 
trajectory parameters. However, this leads to several transcendental 
equations which can not be solved explicitly for the desired quantities. 
As an alternative, the following indirect approach is considered. 

Prom Chapter II, the equation of the warhead trajectory, over a 
spherical, non-rotating earth is 


r " ^(!-e 


hi 


W 


COS 


4J) 


(3-D 
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where 


is the dlstanoe from the oenter of the earth 


h^ is the angular momentum per unit mass 
e^ is the eccentricity of the ellipse 

is the angle measured from the warhead polar axis 
Similarly, the equation of the interceptor trajectory, ooplanar with 
the warhead trajeotory, is 

hi _ 

^ A * \ 

(3-2) 


r — 


^0-e m ccs 4> m ) 


Now at the point of intercept, the radii of "both the intercept and 
warhead trajectories must he equal, or 






<NX 


= r i 


(3-3) 


where r^ is defined as the distance from the center of the earth to 
the intercept point. Furthermore, if we specify that the two tra¬ 
jectories are tangent at the intercept point, the path angles must be 
the same for both, or 


r mI « 


- 




(3-4) 


where fj ts the angle measured from the local horizontal to the tangent 
to the trajectories (positive counterclockwise). This angle, fj* can 
be found by using Eq (2-31). Therefore 


or 


Iran 


r z = 


-e. sin <t>„ x 
cos 4> w:z 


tan 


Xl - 


- fw sin 4>~x 

l-e m cos 4 mJ . 


(3-5) 


(3-5a) 


Rewriting Eq (3-2) for the intercept point, and by using Eq (3-3) 


we get 


r x 


_ h 


2 

m _ 

cos 



(3-6) 
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We can also write Eq (3-2) for the initial conditions, or launoh 
site on the earth's surface, for which oase 

r __ hi 

° t~e m cos <fi me ) (3-7) 

where is the semi-range angle of the intercept trajeotory. 

For a fixed warhead trajectory and intercept point we know e^, h^ 
<£ WI » Tj and Yj* Equations (3-5*, 6 and 7) oontain four variables, e^ 
h , d> ^ and ep . Therefore, if we fix 4> , we have three 

equations and three unknowns and the theoretical intercept trajeotory. 
parameters are completely determined. To find these parameters, we 
first divide Eq (3-6) by Eq (3-7) to give 

r x __ / - e m cos 4> mo 

r ° /_e « cos 4> mx ( 3 _8) 

Now solving Eq (3-5a) for e , we get 

m 

e » -tan Xr __ 

m SLn ^rnx - tan T z cos 4> mX (2,-9) 

Substituting Eq (3-9) into Eq (3-8) and simplifying, 



sin 4> mT - tan 7T X cos 4> mx + tan 4 mX cos 4> mo 

SLn 


(3-10) 


By rearranging Eq (3-10) we can get 


sin 4>mx 
tan r x 


(-* ' ') + 005 


= 003 4>mc 


(3-11) 


Next, let us define a new quantity which is a constant for the 
intercept point. Let 


/ 

tan r z 



(3-12) 
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Substitute Eq (3-12) into Eq (3-ll) and rearrange. 

cos 4> mx + 'S sin 4> mx = cos 4> mo (3-13) 


Equation (3-13) contains only one unknown, <f> , and this can be 

solved as follows? Construct a right triangle as shown in Pig. 3-3. 
By referring to Pig. 3-3 it can be seen that 


/ = sin A (3-14) 

cos A (3-15) 



If we substitute Eqs (3-14) and (3-15) into Eq (3-13) the result is 


'fi+Y r sin A COS 4> mx + yV y* 1 cos A sin 4> mx = cos + mo 


(3-17) 


which can be rewritten as 


By solving Eq 


sin(A + = 

(3-18) for <j> mx , 


cos 4 >„ b 

V7* (3.18) 

and using Eq (3-16), the final result is 
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^mx 


sen 



(3-19) 


With 4 > now determined, we can solve for the intercept trajectory 
parameters e m , and h^, hy using Eqs (3-9) and (3-7)» We have also 
indirectly determined the theoretical launch site of the interceptor 
system since the angle from launch site to intercept point, *£^1’ is 
just (See Pig. 3-2) 


Az ~ +- I 


(3-20) 


Velocity Difference Between Warhead and Interceptor 

With the theoretical intercept trajectory now completely defined, 
we wish to find the difference in velocities between the warhead and 
interceptor, so that we may determine how large a rocket is required 
to perform the velocity-match, i.e. reduce the velocity difference to 
zero. Ignition of the velocity-match rocket causes the actual 
trajectory to depart from the theoretical trajectory, and an approximate 
method of correcting for this departure is discussed later. By referring 
to Eq (2-33) it is seen that the angular momentum per unit mass of the 
warhead or the interceptor can be written as 

Av « r x Kvr css r x (3 ‘ 21) 

and / 

hm = r t \z mx cos r z ( 3 _ 22 ) 

If we now divide Eq (3-21) by Eq (3-22) we get 

fom _ Vmx 

h, ~ v„ (3-23) 
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which can he solved for v to give 

mi 


y — \J _l£! 

** wr h„ 


The velocity difference, A v is given hy 


^ V = ~ 


(3-24) 

(3-25) 


By making use of the fact that v „ and v are co-linear, and 

nx wx 

substituting Eq (3-24) into Eq (3-25) we find the magnitude of the 
velocity difference to be 

Av/= v„0-t) 

It is interesting to note that if the angular momentum of the 
interceptor approaches that of the warhead, A v approaches aero, but 
with the intercept conditions we have imposed, this occurs only when the 
two trajectories are identical ... a highly impractical oase. 



Correction to Theoretical Intercept Trajectory . 

In order to correct the theoretical intercept trajectory for 
the effect of the velocity-match rooket, we must make use of the rocket 
equations of motion which were developed in Chapter II. Equation 
(2-70) can be written in vector notation as 

7 b = V„ -h uln D (3-27) 

Equation (2-76) can also be written in veotor notation as 

- ■?. + v. % + S' (3-28) 

where s' is now defined by 





(3-29) 
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If we let our coordinate system coincide with the point where the 
velocity-match rocket is ignited, then b q » 0 and Eq (3-28) becomes 

= V c t b +S' (3-30) 

Now if we re-write Eqs (3-27) and (3-30) and assume that the motion 
takes place in a gravitational field, where g is constant in magnitude 
and direction, these equations are modified to 

V k =V 0 -buln D ( 3 - 31 ) 


and 


5. = V.t k ~hs' +T5t 


(3-32) 


Let us now refer to Fig. 3-4 and assume that our interceptor 
rocket is in free-flight on the theoretical interceptor trajectory. 

If we did not fire the velocity-match rocket, we would coast to the 
intercept point, IP, and arrive there with velocity equal to v ni . 

But we want to match the warhead's velocity, v^ r , at the intercept 
point. Therefore the velooity-oatch rocket must he ignited some 
time before we get to point IP. Let us assume that the velocity-match 
rocket has a burning time equal to the free-flight time from point A 
to IP. Then, when the rocket arrives at point A with velocity v q , 
the velocity-match rocket is fired so that the exhaust velocity vector, 
u, is parallel to the trajeotory-tangent at point IP. If the effective 
mass ratio, D, of the velocity-match rocket is of the right magnitude, 
then the burnout velocity, v^, will be approximately equal to v , as 
seen in Fig. 3-4. However, Fig. 3-4 also shows that the position at 
burnout is point B rather than that desired (IP). 

Suppose we correct this discrepancy by the method shown in 
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Fig. 3-4 

Vector Diagram of Velocities and Displacement 
■for Theoretical Intercept Trajectory 
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Fig* 3-5* First we extend the tangent line ’’up-range*' to an 
aiming point, AP, so that AP-IP is equal to IP-B. Next we adjust 
the trajectory, so that a free-flight interceptor will arrive 
tangent to point AP with the same velocity, v ffiX , that existed at 
point IP of the theoretical trajectory. Now , if the velocity-match 
rocket is fired at some point A’, where A* is chosen so that the 
free-flight time from A' to AP is equal to t^, the burnout position 
will be at IP. The actual boosted trajectory will be along the 
dotted line in Pig. 3-5* 

If we assume that the burning time is short, we can say that 

ffl «lr x \ (3-33) 


Then the change in direction and magnitude of gravity between points 
A' and IP will be extremely small, and by referring to the velocity 
vector diagram in Pig. 3-5 we can say that 


Vo + 9 ** Kwx 


(3-34) 


uln D - v\ - gt k - v c 
uln D=v,-V mz 


(3-35) 

(3-36) 


However, we want v^ to be equal to v , and since u, and 

v are all parallel we can write Eq (3-36) as 
mi 


uln D = v wz - v mI « A v 


(3-37) 


where 4v is given by Eq (3-26). By solving Eq (3-37) for D we 


find that 


D = e 


Av /u 


(3-38) 
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Fi cjf. 3-5 

Vector Diagram oi Velocities and Displacement 
•for Corrected Intercept Trajectory 
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In Chapter IV it will he shown that if we speoify the thrust to 
initial weight ratio, 0, of a rocket, we can find s' such that 

s'= Hu, P, D) 

With D and s' known we can calculate the corrected intercept trajectory 
parameters. 

By referring to Fig. 3-6 it is seen that 

4 - F=90°-'T X (3-39) 

Note that Yj is negative for intercept points down-range from the 
warhead apogee. We can also show that 


tan 6 — 


s'sin F 

n T - s'cos P 


(3-40) 


or using Eq. (3-39) 


, r >5 COS 

tan o - — — -r ~.—— 
r x - s sin ir x 


(3-41) 


Then with 6 known we can show that the flight-path angle at 


point AP is 

-Ji -s 


(3-42) 


where 6 is always a positive angle and the correct signs must 
he used for yj a and yj* By using the law of sines we find that 


r - £lf LH H 
■ r ‘ , sin S 


(3-43) 


or 


_ _ s ' cos Ik 

~ sin 6 


(3-44) 


With Yt i r T » and v _ known at the point AP we can prooede to 
'la la ml 

find the trajectory parameters. 

Re-writing Eq (3-22) for the corrected trajectory at the point 
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The angle from the corrected, interceptor launch site to the aotual 
intercept point is 

= + ftnxa + & (3-51) 

where 6 is given hy Eq (3—41). 
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In this chapter, first we derived the theoretical intercept 
trajectory equations hy fixing the warhead trajectory and the 
intercept point, and imposing the condition that at the intercept 
point the warhead and intercept trajectories had the same radius, 
r^, and the same path-angle, y^. By assuming a short burning 
time for the velocity-match rocket, and that g was constant in 
the vicinity of the intercept point, we then corrected the 
theoretical intercept trajectory for the effect of the rocket dynamics, 
and obtained the parameters for the corrected intercept trajectory. 

We also determined the effective mass ratio required for the 
velocity-match rocket. With this trajectory completely determined, 
and the effective mass ratio of the velocity-match rocket known, we 
can now consider the boost-phase of the intercept. 
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IV. “ Boost-Phase Analysis 

In Chapter III it was shown that if the warhead trajectory, 
intercept point, and semi-range angle for the interceptor were 
fixed, a unique, free-flight trajectory is determined. We also 
found the required mass ratio for the velocity-match rocket. In 
effect, we are now at the point where we have a "payload*', and a 
desired trajectory, and must find the "booster rocket requirements 
needed to get this payload onto this trajectory. 

There are an infinite number of points where we can enter this 
intercept trajectory, and for each point the required burnout 
velocity and direction is different. This point must also be at 
some altitude above the earth because of the time required to 
accelerate to the final desired boost-phase burnout velocity. If 
we also consider the fact that most rocket launchings consist of a vertical 
flight segment, a segment at constant angle of attack, a segment at 
zero lift, and a constant attitude segment, it can be seen that there 
are many possible flight paths from the launch site to the burnout 
point. The problem is further complicated by atmospheric drag, 
variable gravity, and the variation of rocket thrust with altitude. 

Thus there does not appear to be a neat, closed-form solution to the 
problem, although many approximate methods have been devised (Ref 12). 

In WADC TR 57-724 (Ref 2) a method for the preliminary 
design and optimization of a long-range, ballistic missile is 
outlined. Basically, it is a trial and error method, and while it 
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is well suited to the design of one particular missile, it does 
not lend itself to this investigation, shore we have many 
trajectories to consider. 

Approilmations 

To proceed with the problem, let us make some assumptions 
which enable us to develop a method for approximating the required 
thrusts and weights, and the available warning time for the 
interceptor system. First we assume that the boost-phaee is a small 
portion of the intercept trajectory, allowing us to consider gravity 
as a constant vector. This enables us to choose an inertial, 
rectangular coordinate system on the earth’s surface near the 
launch site, as in Fig. 4-1* Next we neglect the vertical, tip- 
over, and constant angle of attack segments of the boost-phase, 
and assume that the thrust vectors of each stage are applied at a 
constant direction relative to the inertial coordinate system. 
Although this is not practical for an actual launching, this is a 
reasonable approximation since the vertical and tip-over segments 
are usually a small portion of the boost-phase. Then we assume that 
the aerodynamic drag is small compared to the booster rocket thrust, 
and can be neglected. Finally, we neglect the variation of thrusts 
due to altitude, and assume that the specific impulse and structural 
factor are the same for each stage. 

The K-Stage Boost Socket 

With these assumptions in mind, let us develop seme equations 
for the k-stage rocket as shown in Fig. 4-2. Some useful relation- 
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ships for this rocket are listed as follows (Ref 3)* The 
symbols used are shown in Pig. 4-2. 

Structural factors 

*tll 

initial mass of the j stage 
n, B — - ' —— 

burnout mass of the stage 


or 




effective mass ratio of the 


.th 

0 


sub-rocket s 


(4-1) 




_ Mj _ nr\j + m }H j- ■■■■+■ rn k -f- m? _ 

Mjf -^r ■+- m j*' + m s+i + " + "V ^j~(~n~) m j 


(4-2) 


where M . is the burnout mass of the j** 1 sub-rocket. 

Jf 

stage mass ratio: 




(4-j 


overall mass ratio: 


r > 






(4—4a) 


or 


Y.- 1 +% + *' + ■■■+% 


(4-4b) 


Optimization . It can be shown that the optimum arrangement of 
sub—rocket masses, which results in a maximum final burnout velocity 
(vacuum conditions, constant g), for a rocket of constant overall 
mass, with the same exhaust velocity of each stage, is given by the 
equation (Ref 3) 


0, 




(4-5) 
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If the structural factors of each stage are the sane then the equation 
for optimization becomes 


D, = D 2 = • • • ~ D k (4-6) 

The overall mass ratio for optimum staging, with the stage structural 
factors and exhaust velocities equal, is 


K -/W 


(4-7) 


The sub-rocket mass ratios are then 

% = (K - 



(4-8) 


Squations of Motion. In Chapter II it was shown that the 
equations of motion of a one stage rocket in field-free space are 

v k =v 0 + ulnD (2_ 7 o) 

+ r ( 2 - 76 ) 

If we re-write these equations in scalar form for our rectangular 
coordinate system and consider the effect of gravity they become 


x b = k 0 + u cos X In D 

(4-9) 

= + u SLn Xln D -gt b 

(4-10) 

x c t k + s' cos \ 

(4-11) 


y.= y 0 t k + s sin A 

where s =0 
o 

= the horizontal component of burnout velocity 
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y^ ■ the vertical component of burnout velocity 
x q = the horizontal component of initial velocity 
y Q = the vertical component of initial velocity 

X = the angle between the thrust veotor and the horizon 


and 



( 2 - 76 ) 


Let us now find the components of burnout velocity and displacement 
for a k-stage rocket, starting from rest. 


Velocity: 


Stage 1. 

X, = u cos A In D. 

b i 

(4-13) 


y bi = u sin A In D, -gt fc 

(4-14) 

Stage 2 

k h = k h + u cos A In D 2 

(4-15) 


Stage k. 


y b = y ki ~h u sin A In D 2 - gt ^ ( 4 - 16 ) 

X, = X, -h u cos >Zn 0 K (4-17) 



\ = y v , -t- u SLn * ln D * -9% 

(4-18) 

Displacement: 



Stage 1. 

X b = -s' cos A 

(4-19) 


, • x sK 

y h) = s, sm A - — ' 

(4-20) 

Stage 2. 

\=\+\ \+K cos A 

(4-21 i 


- y», +■ y 4 , \ + h st ” * 

(4-22) 
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Stage k. y = x . ~h X, t, -h s' COS A 

k K-< k *-/ k K k 


(4-^3) 


>4 , +->V *V, +- K sinb- 


9 < 


’*-/ " k-l K 

If we now assume that we have a rocket with optimum staging 


(4-24) 


(D's all equal), and further assume that the burning times, t. , 

D 

for all stages are equal, we can say that 


y^ = 2s'sin k y-^t-,-2^ 


Kbw from Bqs (4-27) and (4-28) 

\ + \ + '" + V 
n , + K-i--+K. 

But 


\ln D 

(4-25) 

~2gt k 

(4-26) 

\ In D 

(4-27) 

-kgt, 

(4-^ 


( 4-7 ‘ 


(4-30) 

■ t +"' + Xk )t b 

(4-31) 


(4-32) 


(4-33) 

y h [l + 2-h---+(K-l)j 

(4-34 • 

7 _ k/fc-O 
/ 2 

(4-35) 
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or 






(4-41) 


•j ™j 


By combining Eqs (4-40) and (4-4l) we find that 


t 


= ^_ (3tha 

k J 3 c ft Mj 


(4-42) 


But from Eq (4-2) we can show that 

M, Dj 

Therefore, by combining Eqs (4-43) and (4-42), we find that 


(4-43) 


t 


i; 


“J /Dal) 
Sc & ( Dj J 


(4-44) 


We can also show that 


/ _ Uj / Dj-l-lnDj ) 
J 9 c (3j { Dj / 


(4-4 , >) 


Prom Eq (4-44) it is seen that if the exhaust velocities, sub-rocket 

effective mass ratios, and thrust to initial weight ratios are all the 

same, then the burning times are also all the same. This is 

exactly the assumption we made to get Eqs (4-27» 28 , 36 and 37). By 

using Eqs (4-44) and (4-45) it can be shown that the final velocity 

th 

and displacements at burnout of the k stage are 

X. = kuln D cos X ( 4 - 46 ) 

6 k 

y k ^ = ku In D sin X - (4-47) 

ku*cosX_ rZCD'l-ln D) + (k-0(P-Qln D 1 
£ S */S j D J (4-48) 
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ku z sin A f 2(D-t-ln D) + (k-0(D-!)ln 0 

2 9c/ £? 7 D J 


/ rku(D-l) T 
2 9c[ P D J ( 4 - 49 ) 

Thus, with the burning time eliminated from the equations hy use 
of the thrust to weight ratio, we oan say that 

Burnout Conditions « f(k, u, X, j5, D) (4-50) 

We must now combine these equations with the correoted trajectory 
equations of Chapter III, so that we may find a method for evaluating 
a function of the form 

I) « f(k, u, X, p, Burnout Conditions) (4—51) 


Technique for Satisfying Burnout Velocity Requirement 

First let us solve the problem hy satisfying only the conditions 
of burnout velocity and direction, and later consider what must be 
done to also satisfy the condition of burnout position. Assume that 
we are given a required burnout velocity, v^, and burnout angle, 
as shown in Fig. 4-1» and that we want to find the required sub-rocket 
mass ratios, D. We assume also that k, and £3 are given. We can say 

that 


tan 77 
* y 



(4-52; 


where and y^ are given by Eqns (4-46) and (4-47) with the 
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k subscripts omitted for simplicity. Then 

tan Tj t = tan A - 

If we multiply both sides of Eq (4-53) by the quantity, cos A oos r) 
and rearrange, the resulting equation is 

\ \ (D-l) cos il* 

stn A cos 7f k - cos A sin 7^ - —^l ^ D — (4-54) 


Then 

A \ (D-l) cos 7? a 

sin(A- Vk 

(4-54a) 

Solving Eq (4-54a) for A, 



V • -ir(D-l) cos TUI . 

sm L fiD VTd~J + X 

(4-^>) 

or 

\ - ??, +■ sin'' f 

(4-5:) 

where 

K> (D-0 COS Tl k 

S fiD In D 

(4-57) 


From Fig. 4-1 the horizontal component of burnout velocity is 

h = V h cos V k (4-58) 

By combining Eqs (4-46, 55 > 57, and 58) we can show that 

ku In D cos (y k -h sin' 1 f) - v k cos nl h =■ O (4-59) 


or -F(D)= O 

Now we have an equation with one unknown, D, which would be 


(4-60) 
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difficult, if not impossible, to solve explicitly for D. A numerical 
solution is found by using Newton's iteration method (Ref 14*495)• 

With this method, we choose a trial value for E, say For this 

value we evaluate f(B^) and f'^) where 


The equation for f'(B) can be shown to be 

rrD)*ku^ 


(D-l) cos A f (V7 ~ In D) sin A 

D(D-1) — 


(4-6la) 


Then our new trial value for H is given by the formula 



£ (Dj) 
■f'(Oi) 


( 4 - 62 ) 


This process converges quite rapidly in this particular case, and 
is repeated until the required degree of accuracy is attained. Thus 
we can find D, the sub-rocket effective mass ratio required to 
achieve a given burnout velocity and direction. 


Technique for Satisfying Burnout Position Requirement 

In order to simultaneously satisfy the oondition of burnout 

altitude, y^, or radius, r^, we must resort to a second trial and 

error process. To do this, first choose a trial value for burnout 

range angle, 9^, (See Fig. 4-3a). With the chosen 4> ^ and the 

corrected interceptor trajectory parameters, we use the general 

trajectory equations of Chapter II to determine r v , v, , and v, . 

°1 °1 °1 

As a first approximation let us assume that the launch site in the 
rectangular coordinate system coincides with the theoretical launch 
site in the polar coordinate system. Then we can say that 
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^nioa 


(4-63) 


where.these angles are shown in Fig, 4-3a. Therefore the first 
trial hurnout angle in rectangular coordinates is 

^ = %-€, (4 " 64) 

Then by using Eq (4-59) we can find a value for With known, 

Eqs (4-48) and (4-49) are used to find x, and y, . 

°1 b l 

We can now compute the actual value of r^, which is designated r^ . 
using the relationship (See Fig. 4-3a) 



= V^+yJ 


* + 4* 


(4-65) 


Now we can remove the restriction that the rectangular and polar 
launch sites are coincident by finding that 


= tan 


-I 


r n +■ 


A 


( 4 - 66 ) 


and change from the coordinate system of Fig. 4-3a to the adjusted 

coordinate system of Fig. 4-3b. With the value of r. from Eq (4-65) 

°2 

and the interceptor trajectory equations we can find new values of 
4 >, , v, . and v. . Now substituting €. and v, into Eq (4-64) we 

b 2 V *2 2 *2 ’ 

find ti^ . Then again using Eqs (4-59? 48, and 49) we find values 

for Dg, and y^ . The iteration process is continued until 




(4-67) 


Note that we use the previous value of 6 to compute a new value 
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for rj . This introduces a small error for the first few trials 
b 

since in all likelihood the relative positions of the rectangular 

and, polar launch sites will be different from one trial to the 

next. However, since this process is convergent, as the number 

of trials increases, € will approach a constant value and the 

error will approach zero. An exaggerated illustration of how this 

iteration process might converge is shown in Pig. 4-4* Thus, by- 

using this double iteration process, we find the effective sub-rocket 

mass ratio, D, for the boost-phase rockets, and values for 4> and r, . 

b b 

We are now prepared to achieve our objective, and determine 
the weights and thrusts required, and calculate the excess time. 

First let us find the weights and thrusts. 


Determination of Weight and Thrust Requirements 

The effective sub-rocket mass ratio for the velocity-match 
rocket is found by using Eq (3-38). Let us also specify a thrust 
to initial weight ratio for this rocket, and call this quantity 
P u , where the u subscript is used to avoid confusion with the 
boost-rocket parameters. Assuming a one-stage velocity-match 
rocket, the overall mass ratio of this rocket is found from Eq (4-7) 
and Pig. 4-2 to be 



m' p _ (n-Q Du 

(n-D u ) 


( 4 - 68 ) 


where ^ is the non-dimensional mass (or weight) ratio of the 
velocity-match rocket (lb of int. rocket per lb of int. vehicle). 
The thrust required for the velocity-match rocket is found from 
Eq (4-38) to be 


= ft 


(4-69) 
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or in non-dimensional form 

Th u - p u % (4-70) 

where T y - the thrust of the velocity-match rocket 

- the weight (or lb mass) of the velocity-match rocket 
Th^ « the lb of thrust required per lb of interceptor vehicle. 
The mass ratios of the boost-stages are found by using Eqs (4-7) 
and (4-8). These mass ratios are 




(4-7) 


(4-71) 

- ( r ° wf 1 vir. J 

(4-72) 


<.4-73) 


The thrusts required for each stage are, from Eq (4-38) and Fig. 4-2 


T, = (3 M, = p(m, + m 2 +'- + m k +m' f ) 

(4-74) 

@(m 2 -h m 3 -f — + m k j- 

( 4 - 75 ) 

T k = ?M k = p(m K + m,) 

( 4 - 76 ) 


By dividing through by m\ and using Eq (4-68) we can write 
these equations in non-dimensional form. 
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Th, - & r„(r, + % + ■■■ + r K + 1 ) 

(4-77) 

Th, = / sr u (%-h% + -+r K + i) 

(4-78) 

77 = + 

(4-79) 


By using Eqs (4—4b) and (4-8) we can further reduce these equations 
and show that 


Th, = 

(4-80) 

-P, _ 

(4-81) 

T ^=/f k t j 

(4-82) 




Therefore, the mass (or weight) and thrust ratios are determined, 
and numerical values of weights and thrusts can he found if m^, 
the weight of the interceptor vehicle, and k, the number of 
boost-phase stages are known. 

Excess Time 

To find the excess time it is necessary to determine the flight 
times to the intercept point for both the warhead and the interceptor 
The difference between these times is defined as the excess time. 

Let us assume that the burning time of the enemy booster 
rockets can be neglected and that the total flight time of the 
warhead from its launch site to the intercept point is the same as 
its free-flight time. This time can be found by using the warhead 
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trajectory parameters and evaluating Eqs (2-46) and (2-47) at the 

* 

angles eh » and <h as shown in Fig. 4-5* Let this time he T . 

The time of the interceptor from its launch site, LS^, to the 

boost-phase burnout point, BO, can be found by multiplying Eq (4-44) 

* 

by K, the number of booster stages. Let this time be T^ . From 

Fig. 3-5 it can be seen that the flight-time of the interceptor 

from point A' to the intercept point, IP, is the same as the free- 

flight time (without velocity-match rocket ignited) from point A* 

to the aiming point, AP. Thus the free-flight time from point 

BO to IP (Fig. 4-5) can be found by using the intercept trajectory 

parameters and evaluating Eqs (2-46) and (2-47) at the angles ep , 

* * 

and 4> . Let this time be T . Then the excess time, T , 

r inxa c x 

is given by the equation 


* , * * x 
T - (T, + T ) 

w b c 


(4-83) 


Therefore, by assuming that the enemy missile is detected at the instant 

* 

it is launched, the interceptor must be launched T^ seconds later, if the 
warhead and interceptor vehicle are to arrive at the intercept point 
simultaneously. 


Maximum Acceleration 

An important consideration, especially if a manned vehicle is used 
in the interceptor system, is the maximum acceleration imposed upon the 
system. There is a peak acceleration at the instant of burnout of 
each sub-rocket, and this is found by using Newton's second law to be 
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where 


v . 
J 


T. 

J 


the maximum acceleration in sea level "g's" 

tb. 

the thrust of the j stage 

"til 

the burnout weight of the j stage. 


But. from Eq (4-38) 


Tj - /Sy Mj 


(4-85) 


By combining Eqs (4-84) and (4-85) it is found that 

4 - Mj 






But 


My 


M 


'--Dy 


By combining- Eq (4-2) ar.d (4-86) we get 

ft = ft ft 


( 4 - 86 ) 


(4-2) 


(4-87) 


Since we have specified that- the ?’s and P’s of each sub-rocket 
are equal, Eq (4—87) shows us that the peel: accelerations during 
the boost-phase are the same. Similarly, the maximum acceleration 
during the velocity-match phase would then be 


\>U - D u 

During an actual launching there would be other components of 
acceleration (eg. angular and centripetal accelerations during the 
tip-over maneuver) but these are considered to be of secondary 
importance for this investigation. 



Summary 

In this chapter a method for determining the booster-rocket 
requirements was developed. This was accomplished by making several 
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simplifying assumptions, developing the equations of motion of a 
k-stage rocket and introducing the concept of thrust to weight 
ratio. The equations derived in this manner were then combined 
with the interceptor trajectory equations of Chapter III to 
determine thrust and weight ratios, available warning time, and 
maximum accelerations of the boost-phase and velocity-match 
rockets. 

Because of the many equations, and two iteration processes 
involved, numerical results are very difficult to obtain without 
the aid of an electronic computer. For this reason, the equations 
are programmed on the IBM 1620 digital computer to permit the 
investigation of many different intercept conditions. The 
computer programs used in the analysis are described in Appendix A. 
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V. Procedure for Obtaining Numerical Results 

Primary Considerations 

The overall vehicle weight, the thrust of the first boost stage, 
the excess time, and the maximum "g" forces on the interceptor vehicle* 
are the principal factors to he determined. 

The excess time is certainly the most important factor in the 
analysis, for without adequate excess time, the intercept technique 
is impractical or even impossible. Therefore we wish to determine 
what intercept conditions allow the most excess time, and the effects 
of the boost parameters on excess time. 

The overall vehicle weight is also quite important, since as 
weight increases, the thrust xoquiredto launch the vehicle increases. 
Since there are limits on the amount of thrust available to perform 
the intercept, thrust rather than weight is the factor which we atteiiq 
to minimize in the analysis. 

Finally, the accelerations which exist during the boost phases 
are also an important factor if the system is manned. Man is 
capable of withstanding approximately ten transverse "g’s" or 322 
feet per second per second applied with the pilot perpendicular to 
tne acceleration vector (Ref. 4:14). Therefore we shall indicate 
the intercept conditions which result in accelerations in excess of 
ten *'g’s'*. 

R*3ic Analysis 

The basic analysis consists of fixing the warhead trajectory, 
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intercept point, and the boost parameters and then varying the 
position of the interceptor launch site. This procedure is 
repeated for different intercept points on the warhead trajectory. 
Finally the complete procedure is repeated for other warhead 
trajectories. Once the basic analysis is complete, the effects of 
varying I with n fixed, varying n with I fixed, and varying I 

S 8 S 

and n in combination are determined. 
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It was shown in Chapter II, that for any given range there 
“is an almost unlimited number of trajectories, depending upon the 
amount of total energy available. Thus we must stipulate more than, 
just the range to define the exact trajectory. In order to more 
completely consider the possible intercept situations, three different 
trajectories are chosen for each range. One is the minimum-energy 
trajectory, and the others are a high angle (high altitude) and a low 
angle (low altitude) trajectory. These trajectories can be completely 
defined if we now specify the theoretical launch angle, y 0 > of the 
warhead. For the minimum energy trajectories, we can use Eq. (2-57) 
to determine y o * For the 6000 mile trajectory 

tan Yq - C08< fr° * 0.64278759 = 0.36192092 

m 1 + sin 4> o 1.76604431 

and y * 0.34906588 radian. 

°ME 

For the 4000 mile trajectory, 

tan y » °- 8 354S776 , 0.53919520 

°ME 1.54950890 

and Y a * 0.49450997 radian. 

TiE 

To establish the high and low angle warhead trajectories for 
each range, we can fix the total energy available for the trajectory. 
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This than determines the energy ratio of the trajectory, which in tura 
determines the theoretical launch angle. Since the choice of total 
energy is an arbitrary one as long as reasonable values are used, it 
is more convenient to specify the energy ratio. An energy ratio of 
0,9 is chosen for the 6000 mile trajectory, and 0.8 for the 4000 mile 
trajectory. Then for the 6000 mile trajectory 


tan T Q = 


1 

2 


ER ^ 

if m y - 4( 1 - er ) 


tan 0 

Vi tan 


0.9 

t ] 

1 

h 

3 

• 

CXI 

ON 

• 

0 

0.9 ) 

1.1917534 

((1.1917534/ 



( 2 - 58 ) 


0.58389717 and 0.17126302 


or y o = 0.52849701 radian (high angle) 

and "if 0 = O.I6961984 radian (low angle) 


For the 4000 mile trajectory 


1 

2 


0.8 


0.8 


0.65771029 yl 0.65771029 ) 


- 4 ( 1 - 0.8 ) 


« I.02032565 and O.I96OI585 

or if = 0.79544624 radian (high angle) 

o 

and -g- * O.I9356362 radian (low angle) 


69 


GAW Mech 6l-2 


These values of 4> and y which are listed in Table I, 
are necessary inputs to the computer solution of the boost-phase 
approximations (Appendix A), since the computer is programmed to 
determine all other parameters of the warhead trajectory from 
these values. 


Table I 

Warhead Trajectory Semi-Range and Launch Angles 
for 6000 and 4000 Nautical Mile Ranges 


Warhead 
Range 
(n.m.) 

Trajectory 

Angle 

Semi-Range 

Angle 

(radians) 

Launch 

Angle ' 

(radians) 


Low 

0.58177638 

0.19356362 

4000 

Minimum Energy 

It 

0.49450997 


High 

II 

0.79544624 


Low 

0.87266456 

O.I6961984 

6000 

Minimum Energy 

M 

0.34906588 


High 

tl 

0.52849701 


Selection of Intercept Points and Launch Sites 

To facilitate the discussion of the position of intercept 
points and launch sites, we will use the angular position of the 
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points relative to the warhead apogee in an earth-centered, polar 
coordinate system. This is illustrated in Pig. 5-1 • 


Polar Axis 



Fi g. S- I 

Position of Launch Site an<i 
Tntercept Point Relative to 
Warhead Apogee 


An intercept point is defined by <£ and is understood 

to be on the warhead trajectory. The position of the interceptor 

launch site, </> T <, , is determined in the same manner, but its position 
ra 

vs on the earth's surface. Note that t> is positive clockwise 
and the warhead apogee is at - 0. With these definitions 

established we can proceed with the selection of intercept points and launch 
sites. 


The intercept points to be investigated are arbitrarily chosen 
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so that a representat.ivs sample of points on the warhead trajectory 

is cheeked. These points are approximately equally spaced along 

each trajectory between warhead apogee and atmospheric re-entry point. 

The interceptor launch sites are determined indirectly by 

using the theoretical intercept trajectory semi-range angles as 

inputs to the computer program. Then the approximate position 

of the interceptor launch site relative to the warhead apogee is 

S - <b . The exact position, which is 4> - <d> 

Y wr r so ~wx r* 

is slightly further up-range than the approximate position (Fig. 3-6)> 
hut this is not critical since we are merely trying to select a 
representative group of launch site positions. The up-range 
positions of interceptor launch sites are limited by practical 
considerations; e.g. no launch sites in enemy territory. The 
down-range position of the interceptor launch site is also limited. 

As the angular position of the launch site approaches the intercept 
point, the interceptor trajectory becomes almost vertical and the 
trajectory-match technique is not practical. 

With these limitations in mind, the interceptor launch sites 
are arbitrarily chosen to determine enough points to plot curves 
which show the effect of launch site position on the boost requirements. 

Selection of Rocket Parameters 

The five remaining variables to be investigated are: 

(1) k, the number of stages for the boost rocket 

( 2 ) I , the specific impulse of the fuel 

s 

( 3 ) n, the structural factor 
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( 4 ) P, 'the thrust to weight ratio for the boost-phase 

(5) P , the thrust to weight ratio for the velocity-match 

rocket. 

In order to continue the investigation some particular values of 

these variables must be selected. 

Number of Boost-Phase Stages . For this investigation three 

boost stages were selected as a compromise between the decreased 

reliability of a large number of stages and the large thrust and 

weight requirements for a small number of stages. This reduces the 

number of variable boost parameters to four. 

Specific Impulse and Structural Factor . Specific impulse is 

determined by the choice of propellant. Values of I vary from 

s 

160 seconds for some old-type solid propellants, to 364 seconds 
for liquid hydrogen-liquid oxygen (Ref 10: 112, 113, 312, 313)* 
Instead of selecting a particular fuel, a reasonable value of I , 

9 

2S0 seconds, is selected for the basic analysis. Later in the 
investigation the effect of different values of I g for a few 
typical intercept situations is determined. 

A similar approach is used to determine the value of the 
structural factor since it depends upon the density of the fuel 
and the detailed structural design of a rocket. A value of n equal 
to xen is selected for the basic analysis. 

Thrust to Initial 'weight Ratios . The choice of numerical 
values for p and appears to be completely arbitrary. Then to 
determine the values to be used in the analysis, it is necessary to 
fix the value for one at an arbitrary value, and determine the effect 
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of the other on overall weight and thrust requirements. Once thes- 
effects are known, the most promising values are used in the "basic 
analysis. 


Summai 


This chapter outlines the general approach used to obtain the 


desired data from the digital computer solution of the intercept 


analysis. These results are presented in Chapter VI. 


74 


GAW Mech 61-2 


VI. Results and Conclusions 


A numerical integration was performed using computer program 
CM-12 (Appendix A), for one particular intercept trajeotory. In 
this integration, the effects of drag and the vertical, constant 
angle of attack, zero lift, and constant attitude flight segments 
of a rocket launching are considered. This was accomplished to 
verify the validity of the assumptions made in the "boost-phase 
analysis. For the case considered, the actual burnout velocity 
came within three percent of that required. The results of this 
integration are summarized in Appendix B. 

Since tne results are only good approximations to the actual 
’-alues because of the many simplifying assumptions made, the coaputt: 
data is presented in graphical rather than tabular form. This also 
gives a better overall picture of this intercept technique, and 
allows important trends to be interpreted more easily. Angular 
inputs and print-outs of the computer are in radians and therefore 
all angles are plotted on the graphs in radians. 

In order to avoid complicating some of the graphs, a code 
is used to identify the intercept conditions. For example, 6H$45 s 30 
la interpreted as follows: 

6 H 6000 nautical mile, high angle warhead trajectory 

45 intercept point on warhead trajectory is 0.45 radians 

down range from the apogee 

30 3 emi-range angle of the theoretical intercept trajectory 
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is 0.30 radians 

The theoretical semi-range angle is selected for this code rather 
than the corrected value, because this is the angle used a3 the 
computer input to obtain the numerical results which are plotted. 

Results 

Effect of on Th^. Figures 6-1, 2, 3 were obtained for 
two intercept conditions: 6 H:45s30 and 6H:45*60. Thrust to 
interceptor vehicle weight ratio, p , was varied while the other 
rocket parameters were held constant. These values are 

I ■ 250 seconds 

s 

n = 10 

3 » 1.5 

k = 3 

For the case, 6H:45*60, has a negligible effect on Th^ (Fig. 6-1). 

thus little advantage is gained by using a more powerful engine for 

the velocity-match rocket. The top curve shows that a reduction 

in thrust, Th^, of approximately 3*5^ can be obtained by increasing 

from one to three; however, since the curve appears to be 

leveling off, little will be gained by using a more powerful engine. 

A comparison of the curves in Fig. 6-1 shows that p^ has more effect 

for the intercept trajectory with the smaller theoretical semi-rangf 

angle (.3 vs . 6 ). This is due to the larger difference, Av , between 

warhead and interceptor velocities. A larger decreases the 

distance s' (See Eq 4-45 and Fig. 3-6). This also decreases r T an 

xa. 

therefore we can say that 
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lim (corrected int. trajectory) =-- (theoretical int. trajectory) 
(3 U —* ao 

In other words we can say that as r approaches r T , the total energy, 

la J- 

E, of the intercept trajectory approaches a constant value (E of the 
theoretical intercept trajectory), thus explaining why the curves 
in Fig. 6-1 level off. 

Effect of p and |3 on v. The computer data indicated that for 
the values of |3 used, the maximum acceleration during the boost-phase 
ranged from approximately 3*0 to 4*5 "g’s", and thus this acceleration 
is not considered to be a problem. Therefore no graphs are plotted 
for this data. 

Figure 6-2 shows that the maximum acceleration of the velocity- 
match rocket is a linear function of p . This is also apparent 
from E^ (4-38), 

^ u~ @u 

where is a constant fox' a particular intercept point. Since D u 
is a function of Av, and Av increases with decreasing semi-range 
angles, the maximum acceleration can be a limiting factor on a 
choice of for intercept trajectories with small semi-range angles, 
especially for a manned interceptor vehicle. For the two cases 
shown in Fig. 6-2, it can be seen that for a given the 
acceleration is more than doubled, when the serai-range angle is 
halved. 

Effect of on T . The same reasoning that was used to 

* 

analyze the effect of on Th^ is used for the effect of on T^ 
(Fig. 6-3). Here again, the effect of j3^ is more pronounced for 
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the smaller semi-range angle, and both curves will eventually 
approach a maximum value of excess time. An inorease in from 
one to three results in an increase in T of only about 3 % and 2.5 $ 
for the upper and lower curves respectively. Although this small 
percent change might be important in some cases, we can safely say 
that in general P u is not an important factor as far as time is 
concerned. 

Effect of p on Th^. Figures 6-4 and 6-5 are plotted from 
computer data obtained by varying p and holding the other variables 
fixed as follows: 

I - 250 seconds 

s 

n - 10 

P = 3.0 

u 

k « 3 

Figure 6-4 shows that for each set of intercept conditions there is 
an optimum value of p which results in a minimum value of Th^. 

These minima are 1.8 and 1.6 for the upper and lower curves 
respectively. However, since the curves are quite flat in these 
regions, a value of 1.5 was arbitrarily chosen as a representative 
value of the optimum p. The existence of these optimum values can be 
explained by considering Eq (4-44), which shows that the burning 
times of the boost-phase stages are inversely proportional to p, 
and the equation 


Th , = p % % 


(4-80) 
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As p is first increased, the burning time shortens, and the 

velocity loss due to gravity ( gt^ ) decreases, thus less fuel 

is required to reach a given velocity and decreases. Hov/ever, 

a point is reached where p increases faster than decreases. At 

this point the thrust required, Th n , is a minimum. 

* 

Effect of p on T . Figures 6-5 shows that an increase in p 
results in a larger excess time, and that this time approaches some 
limit asymptotically. The same explanation that was used in 
discussing the effect of P^ applies here also. We can say that 

lim (interceptor flight time) = (free flight trajectory time) 
ft —> oo 

which means that for an infinite acceleration we can attain the 

theoretical launch velocity at the earth's surface and "coast" 

along the trajectory from the ground up. An increase in p from 1.2 

to 1.5 results in increases in of approximately 10$ and 28$ for 

the upper and lower curves respectively. This means that a 

considerable advantage is gained by using a high value of p, for the 

intercept trajectories with larger semi-range angles, , where 

* 

T^ might be critical. 

Effect of n on Th^. The data used to plot the graphs of Th^ vs 
n (Fig. 6-6) was obtained by holding the other variables fixed as 
follows: 

I = 250 seconds 

s 

P u = 3.0 

p = 1.5 
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1.0 1.2 1.4 1.6 1.8 2.0 


Thrust to Initial Weight Ratio of 
Boost-Phase Sub-Rockets (3 

Fig. 6-5 

_ * 

The Effect of (3 on T 
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k - 3 

The curves show that advanced structural design techniques (higher n), 
result in a great reduction in Th^, hut that a point of diminishing 
returns will he reached, and the curves approach a limiting' value. 

This limiting value is unattainable since as n approaches infinity, 
the stage weight is equal to the weight of fuel. This assumes that 
no structure is required to carry this fuel. The structural factor 
has no effect whatsoever on excess time. 

Effect of I on Th . The curves in Figs. 6-7 and 6-8 were 

8 i- 

obtained from computer data with the other variables fixed as follows: 

e u « 3.o 

P - 1.5 

n * 10 

K = 3 

Figure 6-7 clearly shows why great emphasis is placed upon discovering 
new and exotic rocket fuels. For example, an increase in I from 200 
to 225 seconds results in decreased thrust requirements of approximately 

5C$, and an increase in I from 225 "to 250 seconds results in a further 

8 

decrease of approximately 40 ‘fo. 

Effect of I on T . Figure 6-8 shows that T is practically 

8 X X 

a linear function of I . This can also be seen from the equation 

s 

f _— yi— 

b 3iP\ D J (4-44) 

or -a _ f D- 1 1 

e ( d J 
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Thus as I g increases, burning time increases which results in a 
greater velocity loss due to gravity, and therefore causes a 
decrease in T . For the intercept condition 6H:45 J 60, an increase 
in I from 200 to 225 seconds decreases T by approximately 2.5$> 

8 X 

but this small change is insignificant when compared to the 50 $ 
saving in Th^. 

Combined Effect of Changes in I and n. For comparison 

purposes, the effect of combinations of optimistic and pessimistic values 

of I and n upon Th, for the 6H:45*60 intercept are determined, 
s 1 

These results are shown in Table II. 


Table II 

Combined Effects of 

Specific Impulse and Structural Factor 
on Overall Thrust Requirements"' 


Specific 

Structural 

Th l 

Impulse (sec.) 

Fact or 


8 

796 

200 

20 

251 


8 

45 

350 

20 

30 


* Note: For 6H:45s60 intercept, (3 = 1.5. = 3.0 
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Effect of <fi LS and <f> on Th^ and T . Figures 6-9, 10, 
m 


and 11 show the effect of launch site position, TC( , on Th, for 

Lb 1 

m 

various intercept points on the 6000 nautical mile warhead trajectories. 

* , 

Figures 6-12, 13 and 14 are the curves of T vs 9 TO which correspond 

x hb m 

to the Th^ curves. On the Th^ graphs, the dotted lines (. 


) 


indicate t-iie approximate position of lines of constant T (isochrones) = 
The points marked with a star ( * ) are those points where the 
maximum acceleration of the velocity-match rocket reached 10 "g's". 

For all points above the star on each curve, is greater than 10 
"g's". 

Many of the Th^ curves appear to approach horizontal and vertical 

assymptotes. As the launch site is moved up-range from the warhead 

apogee, mQ& approaches 4 1 the interceptor and warhead 

trajectories will eventually coincide and determine the horizontal 

limit for Th,. The vertical limit occurs when 4> is equal to zero. 

1 ' mo 

However, this is simply a vertical trajectory. The interceptor 

velocity, v , at the intercept point therefore approaches zero as 

<^> approaches zero, which results in a high Av, and consequently a 

much larger velocity-match rocket is required. Thus a much heavier 

boost-phase rocket with higher thrust is required for these cases. 

Illustrations of the Use of the Graphs . To illustrate the use 

of the curves, let us suppose that we have a 10,000 pound ( m ) 

maimed interceptor vehicle, and desire to intercept a 6000 n.m. 

minimum-energy warhead trajectory at a point 0.45 radians ( wx ) 

down—range from the apogee. We also select a launch site that is 0.23 

radians ( <£ ) up-range. By entering Bug. 6-10 at -0.23 on the 

m 
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Launch Site Position Relative to Warhead Apogee <p LS (radians) 
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Fig. 6-11 

The Effect of < Pls„ on Th^ for Various Intercept Points on 
6000 n.m. High Angle Trajectory 
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Fig. 6-lla 

The Effect of Time, Thrust, and n g n Limitations on Location 
of Possible Launch Sites and Intercept Points 
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abscissa, proceeding vertically upward to the curve labeled 0.45> 
and then horizontally to the ordinate, we find that Th^ = 105- We 
also see that the maximum acceleration is less than 10 g‘s and that 

the excess time is between 300 and 600 seconds. By using the same 

. * 

procedure, and Fig. 6-13» we find that T = 410 seconds. The first 
stage boost thrust required is 



As a second illustration let us suppose that our restrictions 
* 

are - 600 seconds and v * 10 g's, and we want to intercept 

a 6000 n.m. high-angle warhead trajectory. Then by referring to Fig. 
6-11 we see that we are restricted to operate in the area between the 
600 second dotted line and a curve joining the stars. This area 
defines the possible combinations of intercept point and launch site, 
and the corresponding values of Th^ that meet the restrictions. 
Furthermore, suppose we have a 5000 pound interceptor vehicle, and a 
first-stage booster rocket with 1,000,000 pounds of thrust, giving us 
a Th. of 200. We are then further restricted to launch sites and 
intercept points below this horizontal line on Fig. 6-11. Figure 
6-lla shows how these restrictions can limit the choice of intercept 
points and launch sites. 
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By comparing Pigs. 6-9, 10, and 11 it is seen that as we 

progress from the low angle to the minimum energy to the high angle 

warhead trajectories, time and g-limitations become less restrictive 

hut higher values of Th^ are required. The curves of T x (Pigs. 6-12, 

13, and 14) are also useful for determining' the possible launch sites 

and intercept points for a given time restriction. If, for example, 

we detected an enemy missile on a 6000 n.m. minimum energy trajectory 

immediately after it was launched, and we required a minimum of 700 

seconds to launch our interceptor system, then we are restricted to 

use values of 4> and on the portions of the curves above 

~ m 

a horizontal line at 700 seconds (Pig. 6-13). 

Figures 6-15 through 6-20 show basically the same information 

for the 4000 mile warhead trajectories that was shown in Figs. 6-9 

through 6-14 for the 6000 mile trajectories. The dotted-line 

* 

isochrones are for equal to 300 and 450 seconds. The absence of 

stars on some of the curves indicates that 10 "g's" was never exceeded 

for these cases. These curves show three important differences for 

the shorter range warhead trajectory : maximum acceleration is not 

as restrictive; lower thrusts are required; and time is much more 

critical. Figures 6-18, 19 and 20 show that for some intercept 

* 

points near the warhead apogee, T becomes negative. This means that 
the interceptor flight time is greater than the warhead flight time, 
•and that the interceptor must be launched before the warhead for these 
cases. Since this is quite impractical we can say in general that 
for shorter range enemy missiles, the intercept points must be located 
further down range, than for longer range trajectories. It is also 
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reasonable to assume that the interception of an 8000 mile missile 

would: impose more severe g-1imitations; require higher values 

of first-stage thrust; and increase the excess time. 

Summary of Results 

1. The thrust to initial weight ratio of the velocity-match 
rocket, |3^ has a negligible effect on overall thrust requirements. 

2. 0 has an appreciable effect on maximum accelerations for 
trajectories with small semi-range angles. 

3- P u has a negligible effect on excess time. 

4. There is an optimum value of thrust to initial weight ratio of 
the boost sub-rockets, |3, for each intercept condition which 
minimizes overall thrust requirements. 

5- For values of p close to the optimum value, maximum 
accelerations are well below human tolerance limits. 

6. A small increase in 0 causes an appreciable increase in excess 
time, especially for small semi-range angle trajectories. 

7« Small increases in both structural factor and specific impulse 
greatly decrease overall thrust requirements. 

8. Specific impulse ha3 a much larger effect on thrust 
requirements than structural factor does. 

9 . An increase in specific impulse causes a negligibly small 
decrease in excess time. 

10. Structural factor has no effect on excess time. 

11. Excess time decreases as the launch site and intercept 
point are moved up-range. 
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The Effect of 9 LSm on T for Various Intercept Points on 
4000 n.m. Minimua Energy Trajectory 
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12. Thrust and weight requirements increase as the launch site 
is moved down-range and intercept point is moved up range. 

13* Thrust and weight requirements increase as the warhead 
trajectory altitude increases. 

14. Excess time increases as warhead trajectory altitude increases. 

15• Higher thrusts and weights are required to intercept longer 
range trajectories. 

16. Excess time decreases as warhead range decreases. 

17. Maximum acceleration increases as warhead range increases. 

18. Time, thrust, and "g" limitations determine the possible 
intercept points and launch sites for a given warhead trajectory. 

19. The validity of the assumptions made in obtaining the numerical 
results was verified for a typical intercept trajectory. This 
makes the graphs a useful tool for obtaining approximate thrust 
and weight requirements for any size of interceptor vehicle, 
for the intercept conditions considered. 

Conclusions 

The trajectory-match intercept technique appears to be 

theoi-etically feasible: 

1. If first stage rocket engines of sufficient thrust to boost the 
interceptor vehicle onto the intercept trajectory are made 
available. Minimizing the weight of the interceptor vehicle is 
a major factor in satisfying this requirement. 

2. If sufficient warning time is available from ballistic missile 
radar or satellite reconnaissance vehicles. 
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3. If low reaction times are possible for the interceptor system. 

One major requirement to meet this restriction is the use of 
solid or storable-liquid propellants. 

4. If guidance systems are developed which can meet the requirements 
of this intercept technique. 

Therefore it is recommended that further study of thi 3 technique 
be made for other warhead trajectories and for different combinations 
of the boost and velocity-match stages. This study could also be 
extended to the more realistic situation of a non-planar intercept 
over a rotating earth. 

In addition, if the interceptor is to be manned, a study should 
be made of the methods which could be used for recovery of the vehicle. 
Two possibilities are a controlled skip-glide re-entry or a parachute 
recovery. 

Even if all of the above requirements are met for this system, 
the final decision as to whether or not this intercept technique is 
the answer to the ICBM defense problem must be based upon an 
operational analysis, in terms of cost, technical and military 
manpower requirements, logistics, and many other factors. Thus, 
the trajectory-match, ICBM intercept technique appears to be 
theoretically feasible, but a more comprehensive study is necessary to 
determine if it is militarily practical. 
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Appendix A 
Computer Programs 

The computer programs used in this analysis are written in IBM 
1620 Fortran Language (Ref 5) for use on the IBM 1620 digital 
computer. The instructions for the use of these programs are given 
in this appendix. The actual Fortran Language Program follows each 
set of instructions. The paper tapes for these computer programs 
are on file with the Mechanics Department of the Institute of 
Technology, Wright-Patterson Air Force Base, Ohio. 

The values of t q and ^ used in these programs are (Ref 7s430): 

(average) =• 3437*747 nautical miles 

p. - 62,628. (nautical miles)// (second) 2 

Program: CM-3 General Trajectory Evaluation 

This program gives the values for r, v, and y for any desired 
position, , on the Ballistic trajectory. It also gives the time of 
flight between apogee and ^ . 

Operating Instructions . 

1. V/hen the computer prints the number (l), the inputs are 4> q 
and of the trajectory. 

2. V/hen the computer prints the number ( 2 ), the inputs are the 
general angle , which defines the point to be investigated, 
and L , which allows the computer to automatically increase 
<p hy A <p after each computation. This gives the option of 
determining the values for regular intervals on the trajectory. 
However, if only one point is desired, make A greater than 
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cfc , since some value of A must be inserted into the computer 
memory. When <f> becomes greater than 4> , the computer sets <fi 
equal to q , calculates the values, and then returns to the start of 
the program (ready to accept new values of < f 3 Q and y q ). 

3* The results are typed in the following order: , T*, r, y, and v. 

4- To change either or A at any time, turn "Sense Switch 1" on. 
The computer will then print the number (l), and new values can be 
inserted in the memory. 


^ ‘ y t» l \'*Z. 




' r 3 T r, S2 ^ 


U-62628. 

RO-3437.747 
5 M-1 
PRINT,M 
ACCEPT,PO,GA 
IF(PO) 1,2,1 

1 E-1.0/(COS(PO)+(SIN(PO))*(COS(GA))/(SIN(GA))) 
HS-U*R0*(1.0-E*(C0S(P0))) 

H-SQR(HS) 

W-1.0-(E**2.) 

G*(H**3•)/((U**2•)*W) 

P-(H**2.)/U 
Z«(SQR(W))/(1,0-E) 

0*2.0/(SQR(W)) 

6 hV2 
PRINT,H ‘ 

ACCEPT,X,0ELTX 

8 S-(-E)*(SIN(X))/(1.0-E*(C0S(X))) 

Y»Z*(SIN(X/2.))/(C05(X/2.)) 

T«G*(-S+Q*(ATN(Y))) 

R«P/(1 .0-E*(C0S(X))) 

A-ATN(S) 

V-H/(R*(C0S(A))) 

PRINT,X,T,R,A,V 
IF(SENSE SWITCH 1) 6,7 

7 IF(X-PO) 4,5.4 
4 X-X+DELTX 

IF(X-PO) 8,5,3 
3 X-PO 
GO TO 8 
END 
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Program: CM-11 Combined. Boost-Phase and Intercept Analysis 

This program contains the intercept and trajectory geometry, 
rocket equations of motion, and boost-phase approximation equations. 
These equations are combined as outlined in Chapter IV to give the 
complete solution to the intercept problem. 

Operating Instructions . 

1. When the computer prints the number (l), the inputs are n, 


2 . 

3. 

4- 

5- 
6 . 


I and (3 . 

8 U 

When the computer prints the number (2), the inputs are 3 and k. 

When the computer prints the number (3), the inputs are 

and v 
'wo 

When the computer prints the number ( 4 ), the input is 


wx 


When the computer prints the number (5), the input is <£> 


mo 


Now the computer solves the problem. The results (for k ■* 3) 
will be typed in the order listed below: 




moa 

*b 

u 


'moa 

kt. 


4 >- 


4> 


mla 


u 


YLX A 

X 


Yo ^1 ^2 Y 3 

Th Th, Th_ Th- 

u 1 2 3 

7. The computer will pause after each line, and the start button 
must be pushed to continue the computation. 

8. When the computation is complete, the operator has several 
options as to which part of the program the computer will go 
to for the next set of intercept conditions to be investigated. 
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options are 

controlled by the 

sense switches 

as follows; 

Change 

SW 1 

SW 2 

sw 3 

^ mo 

Off 

— 

— 

d> 

^ wx 

On 

Off 

— 

(3 or k 

On 

On 

Off 

or v 

wo 'wo 

On 

On 

On 


Program; CM-12 Runge-Kutta Integration 

This program is designed to numerically integrate the boost-phase 
rocket equations of motion for a three stage rocket as a check on the 
results of the boost-phase approximations. 

Operating Instructions. 


1. when the computer prints the number (l), the inputs are I , P, 

s 


and t, . 
b 

2. When the computer prints the number (2), the inputs are Th^, A^, 
Th 2 , A 2 , Thy and Ay 

3. When computer prints the number (3), input t^ (duration of 
vertical flight segment), a (angle between thrust vector and 
velocity vector during tip-over phase), and T (flight path angle 
where tip-over is complete and a is set equal to zero). 

4. When computer prints number (6), input initial conditions for; 
t (usually zero)} <fi (usually zero); r (usually zero); 

Ls (usually - ^ moa ); r (usually r Q ); and y Q (usually TT^). 

3. Turn Sense Switch 1 on. 

o. When computer prints number ( 4 ), input At, the incremental time 
interval of the integration. 
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7. When computer prints number (5)> turn Sense Switch 1 Off. 

Then input At (again), and h, (the factor which determines 
the print-out interval). The results are printed out for 
time intervals equal to hA t, and are in the following order: 
t, <p , r, and v. 
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Fortran Language Program CM-11 


D*3. 

TD«=1 .OE-12 
TR«1.0 
U-62628, 

RE«3^37.747 
PI E-3.1415927 
M»1 

PRINT,M 
ACCEPT,F.Z.BU 
21 M«2 
PRINT,H 
ACCEPT,B,Q 
20 fV3 
PRINT,H 

ACCEPT,POW.GOW 

SPOW*SIN(POW) 

CPOW-COS(POW) 

SGOW»SIN(GOW) 

CGOW«=COS(GOW) 

AA«SPOW*CGOW/SGOW 

EW»1,0/(CPOW+AA) 

HWS=U*RE*(1.0-EW*CP0W) 

HW=SOR(HWS) 

AB-HWS/U 
19 M«4 
PRINT,M 
ACCEPT,PWI 
4 SXA«SIN(PWI) 
CXA*COS(PW!) 

AC«1.0-EW*CXA 

AG«EW*SXA/AC 

RI-AB/AC 

gi=atn(ag) 

CGI«C0S(GI) 

VW1*HW/(RI*CGI) 

16 

PRINT,M 
ACCEPT,PO 

ZZ*((RI/RE)-1,0)/AG 
CPO-COS(PO) 
ZA«S0R(ZZ*ZZ+1.0) 
CAZ-CPO/ZA 
SAZ*SQR(1.0-CAZ*CAZ) 
AZ—ATN(SAZ/CAZ) 
ZB=ATN(ZZ) 

Pl*AZ-ZB 
SPI-SINfPI ) 

CPI»C0S(PIj 
E«AG/(AG*CPI+SPI) 


no 



GAW Mech C, 1-2 


HS-U*RE*(1.0-E*CP0) 

H-SQR(HS) 

VI-H/(RI*CGI) 

DV1-VW1-V I 
GOU/ (RE*RE) 

EV-Z*GC 
DLU-DVI/EV 
DU-EXP(DLU) 

W0U-(F-1,)*DU/(F-DU) 

WU-WOU-1. 

TBU*EV*(DU-1.)/(GC*BU*DU) 
QQ-EV*EV*(DU-1.-DLU)/(GC*BU*DU) 
AJ«(PIE/2.)-G! 

SAJ-SIN(AJ) 

CAJ-COS(AJ) 

TAK-SAJ/((QQ/RI)-CAJ) 
AK-ATN(TAK) 

DP1-PIE-AK-AJ 
SAK-SIN(AK) 

RIA-RI*SAJ/SAK 
G1A*AK-(PI E/2.) 

CGIA-COS(GI A) 

HA-VI*RIA*CGIA 
HAS«HA*HA 
REC-HAS/'J 

ET-(VI*V!/2.)-(U/RI A) 

EAS-1,+(2.*ET*REC/U) 

EA-SQR(EAS) 

AL-U*RE 

AM-U*RIA 

CPOA-(AL-HAS)/(AL*EA) 

CP IA=(AM-HAS)/(AM*EA) 

SPOA-SOR(1,-CPOA*CPOA) 

$PIA—S0R(1.-CP IA*CPI A) 
POA-ATN(SPOA/CPOA) 

PIA=ATN(SPIA/CPiA) 

AFL-POA-PIA+DPI 
TGO«EA*SPOA/(1,-EA*CP0A) 
GO-ATN(TGO) 

PB«0.9*P0A 
DA-0.1*POA 
SPB-SIN(PB) 

CPB-COS(PB) 

12 AN-1,-EA*CPB 

TGB»EA*SPB/AN 

GB-ATN(TGB) 

CGB-COS(GB) 

GBR-GB-DA 

SB-SIN(GBR) 

CB-COS(GBR) 

RB-REC/AN 

VB«HA/(RB*CGB) 

VX=VB*CB 
8 DD-D-1. 
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DL*LOG(D) 

UK»EV*Q 
UU»EV*UK 
QW-Q-1 . 

DB*D*B 

YY«DD*CB/(DB*DL) 

YZ-SQRO ,-YY*YY) 
AP*ATN(YY/YZ) 

A=AP+GBR 

SA=S!N(A) 

CA=COS(A) 

F D*UK*DL*CA-VX 

1 F(FD*FD-TD) 1 ,1 ,2 

2 DFA«DD*CA+YY*(DD-DL)*SA/YZ 
DFD-=UK*DFA/(D*DD) 
D*=D-(FD/DFD) 

GO TO 8 

I AU*=2 ,*(DD-DL)+QW*DD*DL 
AR«UU/(2„*GC*DB) 

XOAR*CA*Ai2 

YOAR*SA*Ai^-AR*Q*DD*DD/DB 
0A=ATN(XC/(YC+RE)) 
SDA-SIN(DA) 

RBA*=XC/SDA 
iF(RB-RBA) 3,7,9 

3 RC-RBA-RB 
GO TO 10 

9 RC-RB-RBA 
10 IF(RC-TR) 7,7,11 

II RB*RBA 

CPB-(1.-REC/RB)/EA 
SPB*SQR(1,-CPS*CPB) 

PB=ATN (SB/CPB) 

GO TO 12 

7 TB«UK*DD/(GC*DB) 

GOO*GO 
PBB=PB 
PIAA-PIA 
AFL l=*AFL 

PRINT,POA,GOO,PBB,PIAA,AFLI 

PAUSE 

DBB=DB 

DBU=DU*BU 

ALF-A 

PRINT,TBU,TB,DBB,DBU,ALF 
PAUSE 

AW*(F-1.)*D/(F-D) 

WOB*AW*AW 
AS*1.-(1./AW) 

AU*B/AS 
THU=BU*W0U 
I F (Q-2.) 13J3.14 
13 W1B=W0B*AS 
W2B= (l.'0B-W1B)*AS 


ns 
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WO»WOB*WOU 
W1«W1B*W0U * 

W2«W2B*WOU 
PRINT,WU,W0,W1 ,W2 
PAUSE 

thi=wo*b 

TH2-AU*W2 
PRINT,THU, TH1,TH2 
GO TO 15 

14 WOB«WOB*AW 
W1 B*WOB*AS 
W2B»(WOB-W1B)*AS 
W3B»(WOB-W1B-W2B)*AS 

wo«wou*wob 

W1«WOU*W1B 

W2»WOU*W2B 

W3-WOU*W3B 

PRINT,WU,WO,W1,W2,W3 

PAUSE 

TH1»WO*B 

TH2-AU*W2 

TH3»AU*W3 

PRINT,THU,TH1,TH2,TH3 

15 pause 

IF(SENSE SWITCH 1) 1?,16 

17 IF(SENSE SWITCH 2) 18,19 

18 IF(SENSE SWITCH 3) 20,21 

END 



Pi -i 


PI-3.1415927 

EL=LOG(10.) 

U-62628. 

FM-6080.2011 

RE-3437.747 

DIMENSION A(4) t E(4),C(4),TH(3).Q(5).X(5).Y(5),S(3) 
AL-O. 

GC=U/(RE*RE) 

M=1 

PRINT,M 
ACCEPT,Z,B,TB 
EV«Z*GC 
ZB-Z/B 
M=2 

PRINT,M 
DO 3 K-1,3 
3 ACCEPT,TH(K),S(K) 

17 NV=3 
PRINT,M 

ArrCOT T\/ A A AO 

rixy w ui i • i t ^ nr» ^ ni 

DD-O. 

AZ*SQR(0.5) 

A(1)»0.5 
A(2)*1,-AZ 
A(3)-1.+AZ 
A(4)-1./6. 

E(1 )-2. 

E (2 )*=1 . 

E(3)-1. 

E(4)*2. 

C(1)-0.5 
C(2)-A(2) 

C(3)-A(3) 

C(4)«0.5 
DO 2 1-1 .5 
2 0(I)—0* 

PRINT,M 
DO 80 1-1 ,5 
80 ACCEPT,Y(I) 

ACCEPT,GM 
L=11 
K-1 
N-1 
M-4 

PRINT, M 
ACCEPT,H 
J=1 

37 V-SQR(Y(3)*Y(3)+Y(5)*Y(5)*Y(2)*Y(2)) 

GO TO (52,53,54),N 

53 0=0. 

TF=GM-Y(4) 

N=3 
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GO TO 4 
54 TA-TF+Y(4) 

GO TO 70 

52 VV-V*FM 

HT-Y(5)~RE 

HA-HT-86.67 

HB-HT-6.0855 

HC-HT-13.487 

HD-HT-25.658 

HE-HT-28.783 

HF-HT-43.586 

I F(HA) 18,18,19 

18 RH-1./(EXP(EL*(2,624+.108l8048*HT))) 
GO TO 20 

19 RH-1./(EXP(EL*(12.+.01142854*HA))) 

20 IF(HB) 21 ,22.22 

21 AM-.18362-.0040095*HT 
GO TO 40 

22 IF(HC) 23,23,24 
23 AM-.15922 
GO TO 40 

24 IF(HD) 25,26,26 

25 AM-.15922+ # 00l8602*HC 
GO TO 40 

26 IF(HE) 27 ,27,28 

27 AM-,18186 
GO TO 40 

28 IF(HF) 29,30,30 

29 AM-.18186-.0028798*HE 
GO TO 40 

30 AM-.13923 
40 AN-V/AM 

IF(AN-.6) 31,31,32 

31 CD=.25 
GO TO 41 

32 IF(AN-1.6) 33,34,34 

33 CD*AN-.35 
GO TO 41 

34 IF(AN-2.9) 35,36,36 

35 CD-I,73-.3*AN 
GO TO 41 

36 CD—,86 

41 D-.5*CD*RH*VV*VV*S(K) 

4 TA-AL+GM 
70 ST-SIN(TA) 

CT-COS(TA) 

SG-SIN(GM) 

CG-COS(GM) 

X(1)-1. 

GO TO (73.74,75),K 

73 TS-Y(1) 

GO TO 76 

74 TS-Y(1)-TB 
GO TO 76 
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75 TS*Y(1)-2,*TB 

76 X(2)«EV*(CT-D*CG/(TH{K)))/(Y(5)*(ZB-TS))-2.*Y(3)*Y(2)/(Y(5)) 
X(3)«EV*(ST-D*SG/(TH(K)))/(ZB-TS)+Y(5)*Y(2)*Y(2)-U/(Y(5)*Y(5)) 
X(4)«Y(2) 

X(5M(3) 

DO 1 l«1 ,5 

WZ«A(J)*(X(I}-E(J)*Q(I)) 

Y(I)-Y(!)+H*WZ 
1 Q(I)“Q(t)+3.*WZ-C(J)*X(f) 

I F(Y(2)) 71,72,71 

71 TG*Y(3)/(Y(5)*Y(2)) 

GM=ATN(TG) 

72 J-J+1 
IF(J-4) 37,37,5 

5 J-1 

IF(Y(1)-TV) 6,7,7 

7 IF(GM-AP) 8,8,9 

9 AL*=AA 
GO TO 10 

8 AL“0, 

10 IF(Y(1 )-TB) 6,11 ,11 

II IF(Y(1)-2,*TB; 97,98,98 

97 K»2 

IF(L-22) 95 , 6,6 

95 L -22 
GO TO 99 

98 ir(Y(1)-3.*TB) 96,93,93 

96 K-3 

IF(L—33) 94,6,6 
94 L-33 
GO TO 99 
93 K«4 
L-44 

99 PRINT,L 

6 IF(SENSE SWITCH 1) 12,13 

12 IV5 
PRINT,M 
ACCEPT,HG.HH 
H«HG 

HP*Y(1)+H*(HH-l.) 

13 IF(Y(1)-HP) 14.15,15 

15 PRINT,Y(1),Y(4),Y(5),GM,V 
HP«Y(1)+H*HH 

14 IF(K-3) 50,50,16 

50 IF(N-I) 49.49.37 
49 IF(D-DD) 38,39.39 
39 DD-D 

GC to 37 

3P ;F(D/(TH(K))-.0001) 51,37,37 

51 r *~ s 2 
GO TO 37 

16 IF(SENSE SWITCH 3) 3,17 
END 
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^moa 

a 

0.50623 

rad 

d> 

r b 

at 

-.41069 

rad 

^mla 

* 

0,18201 

rad 

&LZA 

a 

0.68583 

rad 

\n 

» 

32.498 

sec 

Kt b 

= 

319.87 

sec 

V 

- 

4.1637 

"g" 

V 

u 

a 

4.9178 

’’g" 

X 

= 

0.73021 

rad 

V 

u 

a 

0.76462 


*0 

a 

72.974 



a 

51.872 



r= 

15.000 



a 

4.3376 


Th 

u 

a 

5.2939 


ft i 

m 

109.46 


Th, 

8 

31.654 


Th 3 

= 

9.1534 



By inputting the corrected, semi-range angle ( < ^ > moa ) and launch 
angle (T moa ) for the intercept trajectory, a print out was obtained 
for points along this trajectory, by using program CM- 3 . 

An interceptor vehicle weight of 10,000 pounds (m^) was chosen, 
and by using the above data the boost-stage thrusts, stage weights* 
and overall system weight were determined to be: 

T 2 = 1,094,600. lb 

T 2 - 316,540. lb 
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T 3 - 91,534. lb 

M Q = 729,740. lb 

m 1 - 518,720. lb 

IDg = 150 , 000 . lb 

= 43,376. lb 


The burning times of the three boost 


stages are equal, and. are 


_b = 106.62 sec 

3 

The total impulse can be calculated by multiplying the thrust 
of each stage by the burning time of each stage. 


Stage 1. 

5-84 

X 

10 4 

ton 

seconds 

Stage 2. 

1.68 

X 

10 4 

ton 

seconds 

Stage 3. 

0.49 

X 

10 4 

ton 

seconds 

Total 

8.01 

X 

10 4 

ton 

seconds 


The equations of motion used are (Ref 2:17) 

r - r<t > 2 = -jq-(Tstn 6 -h L cos r-D r sin <r) - 
rep -/- 2r<fi— sos 0 — L sin J ~D r cos ?r) 


where r 

4 > 

M 

T 


radius from center of earth 

angle from launch site to radius vector 

instantaneous mass of interceptor system 

thrust 
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L - lift (neglected, for this investigation) 

« drag 

(a « earth's gravitational constant 
6 >= angle between thrust vector and local horizontal 

Y ■= angle between velocity vector and local horizontal 

These equations are broken down into four, first-order, simultaneous 
differential equations for the computer program (Ref 9'. 110-120). 

Some estimate of the drag 1 is needed for these equations, since 
the shape and size of the interceptor system is unknown. A curve 
of vs M (mach number), (Ref 2:43)» for a cylindrical rocket with a 
40° apex nose cone, was approximated by straight line segments which 
resulted in the following estimate for C^. 

C d = 0.25 0 < M < 0.6 

C d = 0.25 + (M - 0.6) 0.6 < M < 1.6 

C d - 1.25 - (M - 1.6) 1.6 < M < 2.9 

C, « 0.86 2.9 < M 

a 

The drag is given by the equation 

D r =irC„p v s A 

where e aerodynamic drag, lb 

« coefficient of drag 
d 

p = atmospheric density, slug per ft^ 1 

v * velocity, ft per sec 

2 

A ® cross sectional area, ft 
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Now since the coefficient of drag is a function of mach number, 
and mach number is a function of the speed of sound, we need an 
estimate for the variation in speed of sound with altitude. This 
was obtained by fitting' straight line segments to a curve of a 
versus b, where b is the height above the earth (Ref 1:24). By 
doing this we find that 

a - 0.18362 - 0.0040095 b 0 <b < 6.0855 

a = 0.15922 6 .O 855 <b < I 3.487 


a 

= 0.15922 + 0.0018602 (b - 13.487) 

13-487 < b < 

25.658 

a 

- 0.18186 

25.658 < b < 

28.783 

Si 

= 0.18186 - 0.0028798 (b - 28.783) 

oa ”70-> / V S 

|VJ X U X 

A 1 c 

a 

= 0.13923 

43.586 < b 


where a 

is the speed of sound in nautical miles 

per second, and b is 

the height above the earth in nautical miles. 



In 

a similar manner a curve of log^ p versus b (Ref l:2l) 

is 

approximated by 




log 10 p - - (2.624 + 0.1082 b), 

b < 86.67 



log 10 p « -( 12.0 + 0.01143 b), 

b > 86.67 



The remaining factor in the drag equation is A, the cross 
sectional area of each stage. We must obtain estimates for all three 
stages since burned-out stages are jettisoned. To do this, first we 
assume that the propellant density is 0.06 lb per in . This is an 
average of densities for some commonly used solid propellants 
(Ref 10: 312, 313)» With this density, and the data on stage weights, 
the approximate volume of each stage is found to be: 
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Stage 1. . . 5>000. ft 3 

Stage 2* • • 1 j445* ft 

Stage 3» . . 420. ft 3 

With the volumes known we can now assume some dimensions for the stages 
and estimate the cross sectional areas. These dimensions are as 
follows: 



Selected 

Approx. 

Approx. 


Area-ft^ 

Diam. ft 

Length ft 

Stage 1 

75 

10 

67 

Stage 2 

50 

8 

29 

Stage 3 

40 

7 

11 

The velocity-match rocket 

and interceptor vehicle are, 

of course, 


mounted on top of the three boost stages but it is assumed that their 
cross sectional areas are equal to or less than the third boost-stage 
area. 

The approximate equations for speed of sound, coefficient of 
drag, and atmospheric density are built into computer program CM-12 
(Appendix A), along with the equations of motion. With this program, 
and the calculated rocket parameters, a trial and error method was 
used whereby different values of the launch parameters were selected. 
The values of r^, v^, and were then compared with points on the 
desired intercept trajectory printout, to determine new trial input 
values. After several trials, the following results were obtained: 

Launch parameter inputs 

t * 24 seconds 

v 
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a •> -0.0& radian (4*6°) 

T * 0.8 radian (45*8°) 


where t 


is the duration of the vertical flight segment 

is the angle between the velocity vector and the thrust 

vector 

is the angle where tip-over segment ceases, and the 
zero-lift (thrust along the velocity vector) segment 
begins. That is, when y reaches a value of 0.8 rad, 
a is set equal to zero. 


Computer printouts 



CM-12 (actual) 

CM-3 (desired) 

Approx, 

Error 

r b 

3637.12 n.m 

3639.22 n.m. 

0.065C 

*b 

0.44583 rad 

0.43464 rad 

2.5* 

v b 

3.0322 n.m. per sec 

3.1342 n.m. per sec 

3.3£ 


The approximate error in overall weight (or thrust) due to the 
error in required burnout velocity can be found as follows. It can 
be shown that for a k-stage rocket in field free space, 


r y «- 

H, r n-D,e *« 


Lfo-O,)'** 


where is the overall mass (or weight) ratio corresponding to 


is the sub-rocket effective mass ratio required to reach 
burnout velocity v^, 
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^ and D- correspond to burnout velocity v,_. 

^ 

If we use the preceding computer results, and the above equation, with 

> n ■ 10 




v. 


V, 


u 


1 - - 2,7758 

{3 1*5 

3.1342 n.m. per second (desired burnout velocity) 
3.0322 n.m. per second (actual burnout velocity) 

! g . (250X32.2) ■ 1.325 n.m. per second 


it is found that 


h 

T t 


0.8987 


actual weight 
required weight 


This indicates that there is an approximate error in the overall Bystem 
weight of 11$. Since the first stage thrust is simply 0 times the 
overall weight, this means that 11$ more thrust is required to reach 
the desired burnout velocity, or approximately 1,220,000 lb instead of 
1,094,600 lb. 


It is felt that after several more trials, a more nearly optimum 
flight path could be determined, and that the velocity error could be 
reduced. However, these results indicate that the assumptions used 
in the mathematical developsient of Chapter IV do not invalidate the 
numerical results of Chapter VI. In any case, for heavy interceptor 
vehicles, where thrust is probably much greater than drag, the graphs 
in Chapter VI can he used to determine excellent estimates of the 
thrusts required for different interoept conditions. 
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in CooBictu aagijiwriti^« AiLor I'uuaiiiug hie C O t — i l o tioa mo 
Lieutenant in the USAF Reserve he entered actire duty in October 
1954. Hie military assignment prior to his coming to the Institute 
of Technology was as an Instructor Pilot in the Basic Training 
Program of Air Training Conmand. 



e entered Stevens Institute of Technology, Hoboken, 



New Jersey. He graduated from Stevens in 1953 with the degree of 
Mechanical Engineer and received a commission as Lieutenant in the 
USAF Reserve. After being called to active duty in October 1953 he 
completed pilot training and served as a fighter pilot in Japan. 

His assignments prior to his coming to the Institute of Technology were 
as a Ferry Pilot in the Military Air Transport Service, and Flight 
Service Operations Officer. 



This thesis was typed by Joyce M. Kerr. 
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